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ABSTRACT 


The  primary  purpose  of  this  investigation  was  to  determine  the 
nature  and  characteristics  of  a  class  of  large,  weak  ionospheric  irregu¬ 
larities  which  in  certain  circumstances  significantly  affects  long-dis¬ 
tance  HF  radio-wave  propagation.  A  further  purpose  was  to  determine  the 
relative  practical  importance  of  the  diffraction  produced  by  these  inho¬ 
mogeneities,  in  relation  to  the  influence  of  other  mechanisms  such  as 
Faraday  rotation  and  scattering  by  small  irregularities. 

Measurements  of  quasi-periodic  variations  in  phase-path  length  have 
been  made  on  a  2600  km  one-hop  lower-ray  F-region  path.  Variations  of  8 
to  75  wavelengths  are  regularly  observed  during  winter  daylight  hours. 
They  have  a  quasi  period  of  about  20  minutes.  Correlations  of  measure¬ 
ments  made  at  spaced  receiver  locations  indicate  that  these  variations 
in  phase-path  length  are  produced  by  irregularities  having  a  horizontal 
scale  size  of  about  35  km.  Multiple  frequency  measurements  show  a  verti¬ 
cal  scale  size  of  about  15  km,  with  a  maximum  in  the  lower  F  region.  A 
lower  limit  of  65  m/sec  is  determined  for  the  horizontal  velocity  of  the 
irregularities.  A  downward  vertical  velocity  of  30  to  40  m/sec  is  ob¬ 
served.  The  irregularities  limit  achievable  bandwidth  and  spatial  reso¬ 
lution  for  one-hop  HF  ionospherically  propagated  waves  under  quiet  iono¬ 
spheric  conditions. 

A  model  of  the  irregularities  has  been  devised  and  has  been  used  to 
predict  the  signal-strength  fluctuations  which  would  result  from  propaga¬ 
tion  through  a  medium  described  by  the  model.  Two  classes  of  signal- 
strength  behavior  are  predicted  and  observed  experimentally,  depending 
on  the  relative  magnitude  of  the  phase-path  variations  within  the  limits 
mentioned  above.  Signal-strength  variations  of  a  few  dB,  having  a  20  min 
period,  are  observed  when  the  changes  in  phase-path  length  are  small,  and 
are  attributed  to  focusing  caused  by  the  irregularities.  Larger,  more 
rapid  signal-strength  fluctuations  are  also  observed.  These  are  attribu¬ 
ted  to  multipath  interference  between  several  rays  propagating  through 
the  same  irregularity.  Taken  together,  these  rays  make  up  what  would  be 
identified  as  the  one-hop  lower-ray  F-layer  mode,  which  is  seen  to  be 
less  homogeneous  than  had  previously  been  supposed,  for  propagation  over 
2000-4000  km  distances. 
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I .  INTRODUCTION 


A .  PURPOSE 

The  purpose  of  the  work  described  in  this  report  was  to  determine 
the  large-scale  spatial  properties  of  high-f  reqxiency  (HF)  radio  waves 
which  had  propagated  over  an  undisturbed  (quiet)  one-hop  F-layer  iono¬ 
spheric  path.  These  large-scale  properties  limit  the  maximum  antenna 
size  which  may  be  efficiently  employed  in  HF  systems.  A  further  purpose 
of  the  work  was  to  determine  the  importance  of  ionospheric  irregularities 
relative  to  other  mechanisms  in  producing  spatial  distortions  in  HF  waves. 


B .  BACKGROUND 

Ionospheric  communication  is  made  possible  by  the  refracting  region 
of  free  electrons  which  exists  at  altitudes  of  about  70  to  300  km  above 
the  earth.  Figure  1  is  an  example  of  a  profile  of  electron  density  vs 
height  which  might  be  observed  over  a  fixed  point  on  the  earth.  This 
electron-density  profile  is  typical  of  the  midlatitude  winter  daytime 
ionosphere.  Radio  waves  in  the  range  of  about  one  to  thirty  MHz  are 
bent  sufficiently  that  they  will  return  to  the  earth  at  a  considerable 
distance  from  the  point  of  transmission.  Figure  2  shows  the  paths  of 
rays  at  a  fixed  radio  frequency  for  an  ionosphere  which  is  described 
everywhere  by  the  electron -density  profile  of  Fig.  1. 

At  the  point  of  return,  the  waves  may  be  received,  or  they  may  be 
reflected  from  the  surface  to  make  additional  "hops",  often  reaching 
distances  two,  three  or  four  times  those  for  a  single  hop.  While  tra¬ 
versing  the  ionospheric  portion  of  their  path,  these  waves  are  affected 
by  any  variations  in  electron  concentration  (and  hence  in  refractive 
index)  which  they  may  encounter. 

Variations  (with  time  of  day,  season  of  the  year,  and  location  on 
the  earth)  in  the  angle  at  which  the  sun's  energy  is  incident  upon  the 
ionosphere  cause  spatial  variations  on  a  global  scale  in  observed  pro¬ 
files  of  electron  density  versus  height.  These  spatial  variations  can 
contribute  to  small  frequency  shifts  in  the  received  signals  and  to 
small  effective  "tilts"  in  the  refracting  layer.  Such  frequency  shifts 
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Figure  1.  Ionospheric  electron  density  profile.  [From  TR-112,  SU-SEL- 
69-007,  Stanford  Electronics  Laboratories,  Stanford  Univer¬ 
sity,  Stanford,  California.] 
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and  tilts  do  not  generally  cause  difficulty  in  the  operation  of  high- 
frequency  communication  circuits. 

More  important  degradations  in  the  smoothness  of  received  wave 
fronts  (surfaces  of  constant  phase)  can  be  produced  by  localized  varia¬ 
tions  ir<  the  electron-density  profile  of  the  ionosphere.  These  localized 
variations  (ionospheric  irregularities)  scatter  or  refract  the  incident 
waves.  If  the  irregularities  are  moving,  they  cause  the  amplitude  and 
phase  of  the  received  signals  to  vary  with  time.  Variations  of  this 
sort  limit  the  useful  bandwidths  and  antenna  apertures  that  may  be  em¬ 
ployed  in  HF  communication  systems. 

1 .  Features  of  the  Quiet  Ionosphere 

The  gross  characteristics  of  the  ionosphere,  through  which  ray 
paths  joining  any  two  points  on  the  earth  travel,  may  be  determined  from 
a  sweep-frequency  sounding  over  the  path  (an  "oblique  ionogram",  Fig.  3). 
In  such  a  record  the  time  required  for  a  group  of  waves  to  travel  from 
transmitter  to  receiver  is  measured  over  a  wide  interval  of  radio  frequen 
cies  (Barry  and  Fenwick  '65).  Figure  3a  shows  such  an  ionogram  for  a  mid 
latitude  path  on  a  typical  winter  day.  Some  important  features  are  indi¬ 
cated  in  the  figure.  In  the  frequency  interval  labeled  "single  mode", 
only  one  major  ionospheric  path  exists  between  the  transmitter  and  the 
receiver.  The  ionogram  of  Fig.  3a  was  made  during  a  period  of  undis¬ 
turbed  or  "quiet"  ionospheric  conditions.  Figure  3b  shows  the  effect 
that  an  ionospheric  disturbance  produces  in  an  oblique  ionogram. 

All  of  the  experimental  work  in  this  paper  was  performed  under 
the  following  conditions: 

1.  One-hop,  lower-ray,  F-layer  propagation  only. 

2.  East-west,  midlatitude  path. 

3.  Quiet  ionospheric  conditions. 

2 .  Sources  of  Wavefront  Distortion  in  the  Quiet  Ionosphere 

Information  on  the  spatial  characteristics  of  HF  ionospherical- 
ly  propagated  waves  is  available  in  a  variety  of  forms  and  representing 
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Figure  3  Sweep-frequency  soundings--Bearden,  Ark. 
to  Los  Banos,  Calif. 


a  wide  range  of  ionospheric  conditions.  In  this  section  some  of  the  pre¬ 
existing  relevant  information  and  its  relation  to  the  present  work  are 
discussed . 


a.  Traveling  Ionospheric  Disturbances 

Traveling  ionospheric  disturbances  produce  large  ionospher¬ 
ic  irregularities  having  horizontal  spatial  periods  ranging  from  50  km  to 
over  1000  km.  The  disturbances  occur  at  rates  of  a  few  per  day  for  smal¬ 
ler  disturbances  to  a  few  per  month  for  larger  ones.  Munro  (’50,  ’58), 

Georges  (’67),  Chan  and  Villard  (’62)  and  others  have  studied  such  dis¬ 
turbances.  The  wavefront  distortions  which  result  are  reasonably  well 
understood.  The  low  rate  of  occurrence  of  traveling  disturbances  decrease 
their  influence  on  the  performance  of  I!F  systems.  The  observations  de¬ 
scribed  in  the  present  work  were  made  under  conditions  when  no  large 
traveling  disturbances  were  influencing  propagation  over  the  path  (as 
determined  by  the  absence  of  "kinks"  in  ionograms  over  the  path--see  Fig. 
3). 

b .  Multi-mode  Observations 

The  spatial  characteristics  and  motion  of  the  fading  pat¬ 
tern  produced  by  the  interference  of  rays  propagating  over  different 
ionospheric  paths,  that  is,  by  different  modes,  has  been  studied  by  Ames 
(’64).  The  observed  fading  rates  and  motions  of  the  fading  patterns  are 
adequately  explained  in  terms  of  large-scale  tilts  and  motions  of  the 
ionospheric  layers.  The  observations  of  the  present  study  are  limited 
to  the  frequency  range  where  only  a  single  mode  is  propagating. 

c .  Single-mode  Oblique-incidence  Measurements 

Measurements  during  periods  of  quiet  ionospheric  conditions, 
using  frequencies  for  which  only  a  single  mode  is  propagating  (or  where 
modes  are  resolved  using  pulse  methods),  have  been  made  by  a  number  of 
investigators.  Such  measurements  permit  the  most  direct  interpretation 
in  terms  of  causative  mechanisms  in  the  ionosphere. 

One  category  of  single-mode  oblique-incidence  studies  con¬ 
siders  the  scattering  produced  by  large  numbers  of  small  (less  than  5  km 
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horizontal  extent)  weak  irregularities.  The  statistical  properties  of 
the  resulting  field  are  measured  to  determine  the  statistics  of  the  scat- 
terers.  It  is  assumed  that  the  received  wave  is  composed  of  a  single 
specular  component  plus  a  large  number  of  weak  scattered  waves.  The 
scattered  energy  introduces  amplitude  and  phase  fluctuations  in  the  re¬ 
ceived  signal  which  are  statistically  distributed  in  a  characteristic 
manner  that  defines  the  ratio  of  the  unscattered  (specular)  energy  to 
the  scattered  energy.  This  ratio  (p)  is  often  referred  to  as  the 
coherence  ratio  of  the  received  signal.  Al'pert  ('63)  summarizes  a 
number  of  amplitude-distribution  measurements  of  ionospherically  propa¬ 
gated  waves.  Whale  and  Gardiner  ('66)  derive  relations  between  the 
statistics  of  the  phase  difference  measured  on  two  spaced  antennas, 
and  the  coherence  of  the  received  wave  field.  Coherence  measurements 
have  not  produced  a  clear  description  of  the  influence  of  small  irreg¬ 
ularities  on  ionospheric  communication.  Phase  measurements  usually 
indicate  much  larger  values  for  the  coherence  ratio  than  are  found 
using  amplitude  measurements. 

A  ■•cond  category  of  studies  considers  the  effect  that 
a  smaller  number  of  larger  irregularities  has  on  waves  propagating 
through  the  ionosphere.  The  effects  of  each  individual  irregularity 
are  measured,  frequently  at  widely  separated  locations.  These  studies 
have  mainly  been  concerned  with  irregularities  in  motion,  primarily 
because  the  observation  of  temporal  variations  is  experimentally  easier 
than  the  observation  of  spatial  variations.  The  velocity  and  size  of 
large  irregularities  have  been  determined  by  measuring  the  time  separa¬ 
tion  of  the  effects  of  a  given  irregularity  at  spatially  separated  points. 
The  measurements  of  traveling  ionospheric  disturbances  described  above 
fall  into  this  category.  Also,  Bramley  ('53,  '55)  has  studied  variations 
in  the  angle  of  arrival  of  obliquely  propagated  waves  at  sites  separated 
by  27  km  to  213  km.  He  measures  angular  deviations  in  azimuth  of  about 
a  degree,  having  a  pseudo-period  of  about  20  minutes.  These  deviations 
are  highly  correlated  at  27  km,  but  are  nearly  uncorrelated  at  the  213  km 
spacing.  Although  it  is  difficult  to  determine  the  shape  and  strength  of 
irregularities  from  angle-of -arrival  data  only,  Bramley  attributes  the 
observed  variations  in  angle  of  arrival  to  weak  traveling  irregularities 
having  spatial  periods  of  50  to  200  km. 
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A  third  category  of  measurements  of  the  spatial  properties 
of  ionospherically  propagated  waves  involves  measuring  the  amplitude  and 
phase  of  the  waves  at  a  large  number  of  spatially  separated  points.  (The 
multi-mode  work  of  Ames,  described  above,  used  this  technique.)  The  major 
work  of  this  type  for  single-mode  propagation  under  quiet  ionospheric  con¬ 
ditions  is  that  of  Sweeney  ('70).  He  employed  eight  receiving  elements 
equally  spaced  along  a  line  2.5  km  long.  He  observed  relatively  linear 
wavefronts  having  mean  angle-of -arrival  variations  similar  to  those  ob¬ 
served  by  Bramlry.  In  addition,  Sweeney  observed  amplitude  differences 
of  about  10  dB  and  phase  fluctuations  of  about  ±20°  on  the  individual 
elements.  These  fluctuations  are  more  rapid  (having  periods  of  less  than 
1  min)  than  the  variations  in  the  mean  angle  of  arrival.  The  mechanism 
that  produces  the  more  rapid  fluctuations  is  not  determined. 

C.  APPROACH  IN  PRESENT  STUDY 

This  report  describes  a  series  of  experiments  performed  to  determine 
the  large-scale  (large  compared  with  the  2.5  km  aperture  used  by  Sweeney) 
spatial  properties  of  ionospherically  propagated  waves.  A  new  phase-mea¬ 
suring  technique  is  used  to  study  the  shape  of  the  observed  wavefronts. 
Phase  measurements  are  made  using  a  fixed  receiver  site  plus  an  airborne 
receiving  system.  The  use  of  an  airborne  receiving  system  allows  phase 
changes  to  be  measured  simultaneously  at  two  points  whose  separation  may 
be  changed  rapidly.  The  moving  airborne  receiver  is  sensitive  both  to 
moving  disturbances  and  to  wavefront  distortions  which  are  fixed  in  space, 
that  is,  not  varying  significantly  with  time. 

Data  gathered  by  use  of  this  technique  reveal  the  presence  of  a  new 
class  of  ionospheric  irregularities,  which  appear  to  be  a  continuation 
to  smaller  sizes  (approximately  35  km)  and  weaker  strengths  (approximately 
1  per  cent  change  in  electron  density),  of  the  irregularities  observed  by 
Bramley  and  others.  Simultaneous  measurements  with  several  radio  frequen¬ 
cies  are  used  to  determine  the  vertical  size  (about  15  km)  of  these  new 
irregularities.  A  computer  model  is  then  used  to  predict  the  influence 
of  the  irregularities  on  propagation.  Amplitude  fluctuations  which  re¬ 
sult  from  propagation  through  such  irregularities  are  predicted  from  the 
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model,  and  are  compared  with  the  fluctuations  observed  experimentally. 

The  occurrence  of  these  weak,  medium-sized  irregularities  is  shown  to 
coincide  with  signal-strength  fluctuations.  The  statistical  distribu¬ 
tion  of  amplitude  fluctuations  at  the  receiver  is  measured  and  compared 
with  amplitude-fluctuation  distributions  predicted  from  the  model.  The 
predicted  and  measured  amplitude-fluctuation  distributions  are  in  turn 
compared  to  distributions  predicted  by  other  models  of  fading  mechanisms 
such  as  0-X  interference  and  scattering  by  numerous  small,  weak  irregu¬ 
larities.  This  analysis  can  explain  some  ambiguities  in  previous  studies 
which  attempted  to  expla.in  amplitude  fluctuations  using  a  model  based  on 
the  assumption  of  many  weak  irregularities. 

Determinations  of  the  apparent  vertical  velocity  and  a  lower  limit 
for  the  horizontal  velocity  are  also  drawn  from  the  measurements.  Fi¬ 
nally,  the  performance  of  wide-aperture  antennas  in  the  presence  of  large 
irregularities  is  discussed.  Appendix  A  provides  a  detailed  description 
of  the  experimental  method,  together  with  the  accuracy  requirements  of 
this  method. 


II.  PHASE-PVTH  MEASUREMENTS 


The  experiments  described  in  this  report  were  performed  in  two 
series  of  measurements.  The  objective  of  the  first  was  to  establish 
a  general  understanding  of  the  wave  field  of  energy  which  has  propaga¬ 
ted  over  a  one-hop  lower-ray  F-layer  path  through  the  ionosphere.  The 
results  of  this  first  series  of  measurements  showed  that  the  major  dis¬ 
tortion  of  the  received  wave  fronts  appears  to  be  caused  by  ionospheric 
irregularities  having  a  temporal  variation  with  a  period  of  about  20  min, 
a  horizontal  scale  size  of  approximately  35  km,  and  a  maximum  deviation 
of  electron  density  in  the  range  of  0.6  per  cent  to  6.0  per  cent. 

The  second  series  of  experiments,  performed  about  one  year  after 
the  first  series,  was  concerned  with  determining  the  vertical  scale  size 
(determined  to  be  about  15  km),  and  the  range  of  magnitudes  of  the  vari¬ 
ations  observed.  In  addition,  data  from  the  second  series  of  experiments 
were  used  to  study  the  signal-strength  fluctuations  produced  by  the  ir¬ 
regularities,  and  to  determine  any  observable  motion  of  the  irregulari¬ 
ties  . 


A.  CONDITIONS  OF  THE  MEASUREMENTS 

In  order  to  minimize  the  possibility  of  ambiguities  in  the  experi¬ 
mental  side  of  this  work  and  in  the  interpretation  of  the  data,  the  mea¬ 
surements  were  performed  under  conditions  which  were  expected  to  produce 
relatively  smooth  wave  fronts.  The  path  employed  joins  transmitters  at 
Bearden,  Arkansas  and  receivers  at  (or  near)  Los  Banos,  California  (see 
map,  Fig.  46  in  Appendix  A).  The  use  of  this  midlatitude  path  avoids 
equatorial  and  polar  anomalies.  The  great-circle  distance  from  Bearden 
to  Los  Banos  is  2600  km. 

Measurements  were  made  during  winter  daytime  periods.  At  these 
times,  ionospheric  disturbances  are  infrequent,  and  a  wide  range  of 
frequencies  could  be  used  which  propagated  by  the  one-hop  lower-ray 
F-layer  mode  only.  For  unambiguous  results,  the  measuring  technique 
that  was  used  requires  single-mode  propagation.  Frequencies  as  low 
as  18  MHz  (limited  by  the  appearance  of  the  two-hop  F-layer  mode  or 
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thu  one-hop  E-layer  mode)  and  as  high  as  28  MHz  (limited  by  the  appearance 
of  the  one-hop  F-layer  upper-ray  mode)  were  used.  On  most  days  a  range  of 
between  3  and  10  MHz  was  available,  for  which  propagation  was  by  the  one- 
hop  lower-ray  F-layer  mode  only.  During  winter  nights  and  during  the  sum¬ 
mer,  there  are  usually  no  frequencies  which  propagate  by  only  a  single 
mode  over  the  Bearden-Los  Banos  path. 

A  high-resolution  FMCW  ("chirp")  sounder  was  used  to  monitor  condi¬ 
tions  over  the  path.  This  sounder  provides  oblique  ionograms  having  8 
usee  resolution  in  group  time  delay.  Ionograms  were  made  routinely,  at 
intervals  of  a  half  hour  to  an  hour  and  a  half,  to  determine  the  "single¬ 
mode"  frequency  range  and  to  check  for  the  presence  of  large  ionospheric 
disturbances.  Additional  ionograms  were  made  if  the  phase-path  measure¬ 
ments  (described  below)  indicated  that  any  ionospheric  disturbances  might 
be  present  on  the  path.  Any  traveling  disturbances  (as  studied  by  Munro 
and  others)  produce  noticeable  "kinks"  in  the  high-resolution  ionogram. 

No  data  from  the  phase-path  measurements  were  used  for  periods  when  iono¬ 
spheric  disturbances  could  be  detected  from  the  ionograms. 


B.  PHASE-PATH  MEASUREMENT  METHOD 

The  selection  of  a  measuring  technique  to  be  used  for  this  work  was 
made  on  the  basis  of  the  known  characteristics  of  wave  fronts  propagated 
over  long,  single-mode,  one-hop  F-layer  paths — plus  a  desire  to  determine 
the  nature  of  the  mechanism  producing  wave-front  distortions.  The  use  of 
pulse  measurement  techniques  was  rejected  because  it  was  felt  that  the 
resolution  possible  with  ionospherically  propagated  pulses  (approximately 
2  psec  at  best)  was  not  sufficient  to  provide  any  new  information  about 
small  variations  in  wave-front  shape.  Angle-of -arrival  measurements  were 
considered  impractical  in  view  of  the  desire  to  make  measurements  at  a 
wide  variety  of  receiver  spacings.  (Accurate  angle-of -arrival  measure¬ 
ments  require  very  carefully  selected  receiver  sites.) 

The  use  of  continuous-wave  (CW)  measurements  was  selected  for  the 
following  reasons: 
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1.  Measured  phase  changes  may  be  more  simply  interpreted  in  terms 
of  ref ractive- index  variations  in  the  ionosphere,  than  is  the 
case  for  group-path  measurements. 

2.  The  sensitivity  of  the  measuring  system  to  small  phase  changes 
is  limited  only  by  the  stability  of  the  frequency  standards  em¬ 
ployed  . 

3.  Dispersion — i.e.,  the  propagation  of  different  frequency  compo¬ 
nents  of  a  pulse  over  different  paths — need  not  be  considered. 

4.  The  narrow  bandwidth  required  allows  the  use  of  low  peak  power 
and  reduces  the  system's  sensitivity  to  interference. 

The  major  disadvantage  of  CW  measurements  is  that  frequencies  must 
be  chosen  at  which  only  a  single  mode  is  propagating.  The  path  used  in 
the  present  study  allowed  the  use  of  single-mode  frequencies  during  win¬ 
ter  daylight  hours.  The  methods  of  acquiring  and  analyzing  the  data  are 
described  below. 

Changes  in  the  phase-path  length  (P)  of  rays  propagated  obliquely 
over  a  2600  km  ionospheric  path  were  used  to  study  irregularities  in  the 
received  wave  fronts  and  to  infer  some  of  the  characteristics  of  the 
source  mechanism  in  the  ionosphere.  The  measured  rays  propagated  over 
a  "single-mode",  one-hop,  lower-ray  F-layer  path.  In  the  absence  of  ir¬ 
regularities  in  the  ionosphere,  the  resulting  field  at  the  receiver, 
produced  by  such  "single-mode"  propagation,  consists  of  an  ordinary  and 
an  extraordinary  magnetoionic  ray.  As  a  result  of  the  irregular  bending 
of  rays  by  ionospheric  irregularities,  more  than  one  such  pair  of  rays 
may  arrive  at  the  receiver. 

The  phase-path  length  of  a  ray  propagating  from  transmitter  to  re¬ 
ceiver  through  a  medium  of  varying  refractive  index  is  defined  by  the 
relation 


P  =  f  (ids  ,  (1) 

"'RAY  PATH 

where  p  is  the  refractive  index  of  the  medium  and  ds  is  an  incremen¬ 
tal  distance  along  the  ray  path.  P  is  directly  related  to  the  phase 
cp  of  the  received  wave  by 
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where  \  is  the  wavelength  in  the  medium  and  equals  7^/jj,  in  which 
is  the  free-space  wavelength.  If  changes  in  phase-path  length  (AP)  are 
measured  in  free  space  (p  =  1),  may  be  used  in  relating  Acp  to  AP. 

In  this  context,  cp  is  not  periodic  in  2it  radians.  The  measurement  of 
changes  in  P  as  a  function  of  time  is  made  by  observing  the  changes  in 
cp  as  a  function  of  time — i.e.,  by  counting  the  number  of  cycles  (an  in¬ 
tegral  number  plus  a  fraction)  of  the  received  signal  which  occur  during 
some  sample  interval  (usually  one  second),  and  comparing  this  with  the 
transmitted  frequency.  Such  a  measurement  Is  unambiguous  (unlike  a  sim¬ 
ple  phase  measurement).  That  is,  the  resulting  phase  measurement  is  not 
limited  to  the  range  0  to  2jt  radians. 

The  received  signal  is  converted  to  an  audio  frequency  by  a  phase- 
coherent  receiver,  and  it  is  this  audio  signal  which  is  measured.  The 
use  of  an  audio  frequency  much  higher  than  any  possible  frequency  shift 
in  the  received  signal  avoids  any  potential  plus/minus  ambiguity  .  The 
equipment  and  processing  used  in  these  measurements,  and  the  accuracy 
of  the  measurements,  are  discussed  in  the  appendix. 

If  only  a  single  ray  path  joins  the  transmitter  and  the  receiver, 

the  meaning  of  such  a  measurement  is  clear.  Let  the  received  signal, 

E  ,  be  described  in  the  form 
r 

Er  =  E (t)  sin  j^2nfot  +  cp(t)J  ,  (2) 

where  t  is  time,  f  is  the  transmitted  frequency  (carrier),  E(t) 

o 

is  a  slowly  varying  (compared  to  the  carrier  frequency)  "envelope"  am¬ 
plitude,  and  cp(t)  is  a  slowly  varying  phase  angle  measured  with  re¬ 
spect  to  the  carrier.  The  quantity  cp  will  be  referred  to  as  the 
measured  phase  of  the  received  signal.  Again,  note  that  the  range  of 
9  is  not  restricted  to  2n  radians.  E(t)  is  also  measured  at  the 
sampling  interval.  A  one-second  sampling  interval  has  been  found  to 
be  satisfactory  for  the  relatively  stable  signals  observed  under  the 
measurement  conditions. 
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The  received  signal  E  can  be  represented  as  a  vector  of  length 

r 

E  at  an  angle  cp  relative  to  a  coordinate  system  rotating  at  2 rtf 

o 

rad/sec  (Fig.  4a).  If  truly  only  one  signal  is  present,  this  picture, 
together  with  the  meaning  of  the  measurement,  is  clear.  If  the  received 
signal  is  shifted  in  frequency  from  the  carrier  frequency,  e.g.,  if  a 
doppler  shift  is  caused  by  motions  in  the  ionosphere,  then 

cp(t)  =  2rtfrft  ,  (3) 

where  f  .  (commonly  called  "doppler")  is  the  difference  between  the 
d 

transmitted  and  received  frequencies. 

If  the  received  signal  is  produced  by  a  summation  of  several  rays, 
the  interpretation  of  the  measurement  may  or  may  not  be  clear.  The 
phase  measured  is  that  of  the  vector  sum  of  all  components  present. 

Some  possible  types  of  sums  are  shown  in  Figs.  4b  and  4c  and  are  de¬ 
scribed  below. 

Consider  a  received  vector  composed  of  a  large  signal  vector  (Eg) 
plus  a  small  randomly  phased  vector  (E^).  As  has  been  pointed  out,  the 
measurements  made  in  connection  with  this  work  were  carried  out  under  un¬ 
disturbed  ionospheric  conditions,  using  frequencies  for  which  only  a  sin¬ 
gle  mode  was  propagating.  Measurements  by  Sweeney  (*'0)  and  others  show 
that  under  such  conditions,  phase  changes  are  generally  smooth  for  peri¬ 
ods  of  several  minutes,  but  that  small,  more  rapid  (5-20  sec  period) 
phase  fluctuations  of  about  10  deg  are  observed.  Such  fluctuations  might 
arise  from  noise,  interference,  or  scattering  by  small  irregularities  in 
the  ionosphere.  This  type  of  fluctuation  may  be  represented  by  a  small 
random  vector  E^  added  to  the  main  signal  vector  (Fig.  4b).  The  mea¬ 
surement  of  phase  is  still  unambiguous;  however,  a  phase  error  is  intro¬ 
duced.  As  long  as  the  actual  phase  error  introduced  into  the  vector  sum 
is  not  significant  to  the  data  derived  from  the  measurement,  the  effect 
of  such  perturbations  is  not  important,  and  does  not  modify  the  concept 
of  the  measurement.  Actually,  if  one  received  signal  vector  is  larger 
than  the  sum  of  the  magnitudes  of  all  other  vectors,  the  unambiguous  mea¬ 
surement  is  a  meaningful  representation  of  the  larger  (dominant)  vector 
as  long  as  an  error  in  measured  change  in  P  of  ±7\/ 4  (equivalent  to  a 
change  in  cp  of  ±rt/2  radians)  is  acceptable.  While  an  error  of  rt/2 
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radians  might  invalidate  a  simple  phase  measurement,  errors  of  this 
magnitude  will  be  negligible  in  the  unambiguous  phase-path-length  mea¬ 
surements  . 

Now  consider  a  measured  resultant  vector  made  up  of  two  vectors  of 

similar  magnitude,  E  and  E  ,  whose  lengths  vary  with  time  (Fig. 

bl  b  Z 

4c).  Euch  a  situation  might  represent  two  magnetoionic  components  (see 
Chapter  III);  or  propagation  by  more  than  one  mode;  or  "internal  multi- 
path"*  effects  produced  by  ionospheric  irregularities.  In  the  case  of 
the  ordinary  and  extraordinary  magnetoionic  components,  the  relative  po¬ 
sition  and  strength  of  the  two  component  vectors  will  vary  only  slowly 
with  time.  With  the  exception  of  a  complete  cancellation,  the  measured 
phase  and  frequency  shift  of  the  resultant  vector  will  be  a  reasonable 
estimate  of  the  performance  of  either  of  the  two  rays.  For  rays  of  com¬ 
parable  amplitude  which  do  not  behave  similarly  to  each  other,  we  will 
observe  the  variations  in  length  of  phase  path  P  for  the  stronger  ray, 
within  an  accuracy  of  ±A/4.  Since  the  relative  amplitudes  of  the  rays 
may  vary  with  time,  one  ray  and  then  the  other  may  successively  dominate 
the  measurement.  If  such  relative  changes  in  amplitude  occur  slowly 
enough  (i.e.,  if  many  samples  are  taken  between  changes),  and  if  the 
frequency  shifts  fQ  along  the  ray  paths  are  different  enough  for  the 
two  rays,  we  will  be  able  to  observe  an  abrupt  change  in  the  slope  of 
the  phase-path  record  as  the  change  occurs.  Figure  5  shows  the  type  of 
changes  we  would  expect  under  these  conditions. 

For  larger  numbers  of  components,  a  simple  intuitive  understanding 
cannot  be  continued  unless  one  of  the  above  conditions  can  be  approxima¬ 
ted;  i.e.,  unless  one  ray  or  several  rays  with  very  similar  frequency 
shifts  dominate  the  resultant  vector  for  reasonably  long  intervals  of 
time.  Ultimately,  the  analysis  of  the  amplitude  and  phase  measurements 
will  be  limited  by  sampling  theory  and  measurement  noise;  however,  in 
this  work  the  technique  is  applied  only  under  conditions  when  the  simple 
intuitive  approach  is  applicable. 


*The  term  "internal  multipath"  is  used  to  describe  multiple  rays  which 
propagate  from  a  transmitter  to  a  receiver  via  a  single  ionospheric  mode. 
Internal  multipath  is  present  if  ionospheric  irregularities  produce  suf¬ 
ficient  differential  ray  bending  that  rays  cross  before  reaching  the  re¬ 
ceiver  location. 
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Using  the  measurement  described  above,  it  is  possible,  under  certain 
conditions,  to  determine  changes  in  P  for  the  ray  path  joining  a  trans¬ 
mitter  and  a  receiver,  and  conversely,  it  is  frequently  possible  to  deter¬ 
mine  from  the  measurement  record  whether  the  necessary  conditions  exist 
at  the  time  oi  the  measurement.  Several  implications  of  such  measurements 
will  now  be  discussed. 

Variations  in  the  phase-path  length  of  an  ionospherically  propagated 
ray  can  occur  in  two  ways.  From  Eq.  (1)  it  is  seen  that  variations  in  P 
may  be  caused  both  by  variations  in  p  along  a  fixed  path,  and  by  changes 
in  the  position  or  length  of  the  physical  path.  The  frequency  shift  (dop- 
pler)  which  is  produced  by  varying  the  phase-path  length  of  a  ray  is 


fd  =  ' 


_ o  dP 

c  dt 


where  c  is  the  speed  of  light  and  f  is  the  transmitted  radio  fre- 

o 

quency .  By  measuring  the  slope  of  a  phase-path  record,  we  may  determine 
the  doppler  shift  of  the  ray  that  dominates  the  measurement. 

An  additional  and  perhaps  unexpected  relation  for  rays  propagating 
between  two  fixed  endpoints,  through  a  medium  in  which  the  refractive 
index  varies  with  time,  is  derived  by  Bennett  C'67) 


=  -~  / 


RAY  PATH 


This  relation  states  that  the  instantaneous  frequency  shift  is  a  function 
only  of  the  instantaneous  rate  of  change  of  refractive  index  along  the 
present  ray  path,  and  does  not  include  a  contribution  due  to  the  simul¬ 
taneous  change  in  the  position  or  length  of  the  ray  path  which  may  be 
occurring.  This  fact  will  provide  a  more  direct  interpretation  of  our 
measurements  in  terms  of  the  irregularities  which  cause  them. 


More  precise  statements  ox  Eqs.  (1)  and  (5),  which  account  for  the  an¬ 
isotropy  of  the  ionosphere,  include  an  additional  term,  cos  a,  within 
the  integral,  where  a  is  the  angle  between  the  ray  direction  and  the 
wave  normal  (see  Budden,  1961).  When  considering  small  changes  in  long- 
range  oblique  ray  paths,  the  error  introduced  by  neglecting  this  term  is 
negligible . 
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It  is  also  useful  to  consider  the  technological  constraints  which 
will  influence  the  usefulness  of  the  measuring  technique.  The  measured 
changes  in  phase-path  length  are  not  expected  to  occur  at  a  rapid  rate, 
since  frequency  shifts  greater  than  one  or  two  tenths  Hz  are  rare  under 
the  stated  experimental  conditions.  An  error  in  integrated  phase  aris¬ 
ing  from  the  frequency  sources  used  for  transmitting  and  receiving  will 
be  reflected  in  a  like  error  in  the  measured  phase-path  length.  The 
more  rapidly  phase-path  changes  introduced  by  the  ionosphere  occur,  the 
less  rigid  will  be  the  requirements  on  the  frequency  sources  used.  Fre¬ 


quency  errors  (linear  changes  of  phase)  can  be  tolerated  if  a  linear 


trend  may  be  removed  from  the  resulting  data  without  loss  of  information. 


Instability  in  frequency,  however,  will  degrade  any  data  that  make  use 


of  measurements  of  changes  in  P.  Crystal  frequency  standards  having  a 


stability  of  about  10  will  introduce  errors  in  the  measurement  of  P 


at  a  rate  of  about  1000  m/hr.  Portable  cesium-beam  frequency  standards 


having  a  stability  of  10  will  allow  measurements  with  an  error  of 


about  10  m/hr.  For  the  observed  variations  of  the  present  experiments, 


this  gain  in  accuracy  is  crucial. 


C.  MULTIPLE -LOCATION  MEASUREMENTS 


The  effect  of  refractive-index  irregularities  on  waves  propagating 


through  them  is  principally  determined  by  the  strength  (per  cent  altera¬ 


tion  of  the  refractive  index)  and  the  spatial  extent  of  the  irregulari¬ 


ties  . 


The  first  series  of  measurements  performed  in  this  study  was  primar¬ 


ily  concerned  with  determining  the  strength  and  horizontal  scale  size 


(or  sizes)  ol  the  ionospheric  irregularities  believed  to  be  present  under 


quiet  ionospheric  conditions.  Since  the  size  range  of  possible  irregu¬ 


larities  was  not  known  a  priori,  measurements  were  made  at  as  wide  a 


range  of  spacings  as  possible.  The  work  of  Bramley,  carried  out  under 
quiet  ionospheric  conditions,  indicated  that  the  probable  range  of  sizes 


of  irregularities  to  be  encountered  would  lie  between  30  and  200  km. 
Making  use  of  these  observations,  and  considering  some  practical  experi¬ 


mental  limitations,  a  single  stable-frequency  CW  transmitter  was  used, 
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and  phase-path  measurements  were  made  simultaneously  at  two  receiving 
points  separated  by  distances  which  varied  from  50  to  170  km.  This 
provided  two  ray  paths  in  the  ionosphere  which  were  separated  horizon¬ 
tally  by  25  to  85  km  at  the  maximum  height  of  the  rays.  One  receiving 
system  was  located  at  the  Stanford  University  Wide-Aperture  Research 
Facility  near  Los  Banos,  California,  while  the  other  receiver  was  in¬ 
stalled  aboard  a  DC-3  aircraft.  The  aircraft  was  also  equipped  with 
an  accurate  tracking  system  which,  in  conjunction  with  two  line-of -sight 
stable-frequency  transmissions,  allowed  continuous  (though  not  real-time) 
tracking  of  the  plane's  position.  The  equipment  and  data-processing 
method  employed  are  described  in  detail  in  Appendix  A.  The  resulting 
records  show  the  variation  in  phase-path  lengths  between  the  transmitter 
(located  in  Bearden,  Arkansas,  2600  km  east  of  Los  Banos)  and  each  of  the 
receivers.  The  records  from  the  aircraft  are  corrected  to  compensate  for 
the  variation  of  the  aircraft's  position  from  a  constant  radial  distance 
to  the  transmitter  and  from  a  constant  altitude.  For  most  of  the  experi¬ 
ments,  the  aircraft  flew  a  course  running  from  50  to  170  km  north  of  Los 
Banos.  The  course  was  flown  in  both  directions. 

The  choice  of  this  method  of  observation  results  in  the  ability  to 
explore  wave  fronts  at  various  receiving-antenna  spacings  without  the 
necessity  of  establishing  numerous  field  sites,  and  allows  the  detection 
of  fixed  irregularities.  It  does,  however,  suffer  a  decrease  in  the  num¬ 
ber  of  measurements  that  may  be  made  at  any  given  spacing,  for  the  sake 
of  an  increase  in  the  number  of  spacings  that  may  be  employed.  The  re¬ 
cords  show 


P(t) 


P 

a 


r  =  constant 
6  =  constant 


for  the  fixed  receiver,  and 


P(t,0)  -  P, 
b 


r 


constant 


for  the  aircraft  receiver, 
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where 


P  =  phase-path  length, 


r  =  great-circle  distance  from  the  transmitter  to  the  receiver, 


9  =  bearing  from  transmitter  to  receiver,  and 


P  and  P  =  unknown  constants  at  the  start  of  each  measurement, 
a  b 


The  slopes  of  these  functions  correspond  to  the  observed  frequency  shift 
(doppler)  of  the  signal  and  are  given  by 


—  ,  with  r  and  9  constant,  for  the  fixed  receiver, 
dt  ' 


OP  dP  d@  ,  dG 

^t"  +  o6  *  dt  *  w:*-th  ^  and  r  constant,  for  the  aircraft.  (6) 


The  second  term  in  the  doppler  expression  for  the  aircraft  arises  from 


the  aircraft  velocity. 


Two  sources  of  error  in  interpretation  of  these  data  will  now  be 


mentioned  with  reference  to  their  effect  on  our  observations.  The  first 


is  that  the  ray  path  joining  the  transmitter  and  the  receiver  does  not 


lie  in  a  plane  containing  the  transmitter,  the  receiver,  and  the  center 


of  the  earth,  because  of  the  tilts  in  the  ionosphere  produced  by  the  ir¬ 


regularities  we  are  studying.  Therefore,  the  exact  ray  path  is  not  known, 


and  the  measured  doppler  is  that  corresponding  to  a  point  in  the  iono¬ 


sphere  whose  location  is  not  precisely  known.  The  error  introduced  by 
this  deviation  can  be  estimated  by  making  use  of  angle-of -arrival  mea¬ 


surements  made  by  Sweeney  (’70).  He  used  a  broadside  antenna  array  which 
was  2.5  km  long,  and  which  received  over  the  \rkansas-Los  Banos  path. 


Hu  found  root-mean-square  (rms)  angle-of -arrival  fluctuations  of  about 


1/1  deg,  which  implies  that  the  ionospheric  ray  path  has  an  rms  devia¬ 
tion  from  the  great-circle  plane  containing  the  transmitter  and  receiver 


of  about  7  km  at  its  maximum  height. 
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The  second  source  of  error  is  introduced  by  the  aircraft's  motion 
(the  c)P/d0  term  in  Eq .  6).  The  effect  of  this  error  can  also  be  de¬ 
termined  from  the  known  angle-of -arrival  deviation.  For  a  15  m  wave¬ 
length  (20  MHz),  a  1/4  deg  change  in  angle  of  arrival  will  produce  a 
phase  variation  normal  to  the  nominal  ray  path  of  one  cycle  in  3  km. 

At  the  nominal  aircraft  velocity  of  60  m/sec,  this  will  result  in  a 
frequency  shift  of  0.02  Hz.  As  is  indicated  by  the  measured  data  pre¬ 
sented  in  the  next  section,  these  zero-mean  errors  introduce  only  minor 
errors  in  the  results  for  the  values  of  frequency  shift  measured  in 
these  experiments,  and  for  the  number  of  points  averaged  in  computing 
correlations . 

Measurements  using  the  aircraft  receiver  and  the  Los  Banos  receiver 
were  made  during  April  of  1969.  A  total  of  ten  flights  yielding  satis¬ 
factory  data  were  made.  Frequencies  between  20  and  25  MHz  (wavelengths 
7\  between  12  and  15  m)  were  used  for  these  measurements.  Each  flight 
covered  a  path  between  a  point  170  km  north  of  Los  Banos  and  a  point  50 
km  north  of  Los  Banos,  and  required  between  20  and  40  minutes.  On  any 
given  day,  one  or  two  flights  over  the  path  were  made.  If  two  flights 
were  made  on  one  day,  they  were  made  in  rapid  succession.  (The  data  from 
two  of  the  ten  flights  were  later  discarded  because  of  interference  on 
the  tracking  system  frequencies  in  one  case,  and  evidence  on  an  ionogram 
and  in  the  phase-path  measurements  of  an  ionospheric  disturbance  in  the 
other . ) 

A  typical  pair  of  records  from  the  fixed  and  airborne  receivers  is 
shown  in  Fig.  6.  For  each  pair  of  records,  a  constant  doppler  (constant 
frequency  shift,  indicating  a  linear  time  rate  of  change  of  phase-path 
length)  has  been  subtracted  from  both  records.  This  removes  the  effect 
of  large-scale  ionospheric  layer  motions  which  do  not  influence  the  shape 
of  the  wave  front.  The  most  striking  feature  of  the  records  is  the  quasi- 
periodic  variation  in  phase-path  length  seen  consistently  in  our  observa¬ 
tions.  The  time  between  peaks  of  these  variations  ranges  from  15  to  25 
min.  The  peak-to-peak  phase-path  length  variations  seen  in  these  records 
ranged  from  about  8  7\  (approximately  100  meters)  to  about  50  A  (approxi¬ 
mately  600  meters).  These  length  variations  produced  maximum  frequency 
shifts  ranging  from  about  0.05  Hz  to  0,20  Hz.  Analysis  of  these  data 


(see  Section  D  of  this  chapter)  indicated  the  presence  of  large,  weak 
irregularities . 

The  small  size  of  the  measured  frequency  variations  reveals  some 
of  the  experimental  difficulties  which  may  be  encountered  in  studying 
them.  In  order  to  make  relative  phase-path-length  measurements  at  two 
separated  points,  the  frequency  standard  at  each  point  must  be  suffi¬ 
ciently  accurate  and  stable  that  the  desired  measurement  is  not  masked 
by  variations  in  the  standard.  For  variations  of  the  order  of  10  A  per 
20  min,  at  an  operating  frequency  near  20  MHz,  we  would  like  the  stan¬ 
dard  to  introduce  errors  no  greater  than  10  per  cent  of  the  desired 


measurement.  This  means  one  cycle  per  20  min  at  20  MHz,  or  a  stability 

If  the  two  measurements  are  made  only  a 


(&£/f)  of  about  4  X  10 
short  distance  apart,  a  common  frequency  standard  may  be  used,  making 
possible  relative  phase  measurements  that  are  as  accurate  as  may  be  de¬ 
sired.  However,  if  the  effect  being  investigated  produces  only  an  in¬ 
significant  difference  in  length  of  phase  path  at  points  separated  by 
the  small  distances  which  may  reasonably  use  the  same  frequency  standards, 
independent  standards  of  the  above  accuracy  must  be  used  in  order  to 
enable  observation  to  be  made  at  more  widely  separated  receiving  points. 

It  is  only  recently  that  portable  cesium-beam  standards  have  been  avail¬ 
able  which  meet  this  requirement. 

A  second  difficulty  is  encountered  if  doppler  measurements  are  to 
be  made  using  a  Fourier  analysis  technique  (instead  of  measuring  the 
instantaneous  slope  of  the  phase-path  variations).  The  smallest  fre¬ 
quency  interval  A£  which  may  be  resolved  by  Fourier  analysis  is  1/T, 
where  T  is  the  length  of  the  data  sample.  If  ten  independent  frequency 
intervals  are  desired  in  the  range  of  dopplers  produced  by  weak  irregu¬ 


larities,  i.e.,  f^  =  0.05  Hz,  then  Af  =  0.005  Hz,  and  T  =  200  sec 


During  this  200  sec  interval,  for  the  observed  variations,  the  frequency 
will  have  shifted  through  several  frequency-resolution  cells,  producing 
a  poorly  defined  spectrum.  If  only  a  single  signal  component  is  present, 
we  may  use  a  shorter  analysis  time  and  interpolate  between  frequency  in¬ 
tervals;  however,  this  condition  is  rarely  met  for  ionospherically  propa¬ 
gated  signals.  When  multiple,  comparably  strong  (within  an  order  of  mag¬ 
nitude)  components  are  present,  the  doppler  method  can  separate  them  only 
when  the  frequencies  change  slowly  enough  to  be  resolved. 
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I>.  HORIZONTAL  SCALE  SIZE  OF  THE  OBSERVED  IRREGULARITIES 

An  estimate  of  the  size  (horizontal  extent)  of  the  irregularities 
being  investigated  was  made  from  the  above  data  by  computing  a  correla¬ 
tion  coefficient  ''or  the  variation  of  the  doppler  from  the  mean,  measured 
at  various  spacings  of  the  two  receivers.  The  doppler  was  measured  graph¬ 
ically  from  the  phase-path  records.  A  total  of  eight  pairs  of  records 
was  found  suitable  for  use  in  this  correlation.  For  each  pair  of  records, 


PA  (Ax) 


Ldl 


'd2 


(7) 


was  calculated,  where 


o^(Ax)  is  a  sample  of  the  correlation  for  the  ith  flight  at  a 
receiver  spacing  of  Ax, 

f  is  the  scaled  doppler  value  for  the  fixed  receiver, 
dl 

frf2  is  the  scaled  doppler  value  for  the  airborne  receiver, 
and 

2 

^dl  *’s  mean-scluareci  doppler  for  the  fixed  receiver  for 
the  length  of  the  record . 


The  estimate  of  correlation  was  then  calculated  as 


^  p.(Ax) 

P(^x)  =  2  § 

i=l 


(8) 


The  calculated  points  are  plotted  in  Fig.  7.  The  error  bars  indicate  the 

2 

estimated  variance  (s')-  of  the  correlation  [as  defined  by  Eq.  (7)3 
calculated  from  the  measured  values  by 


(s') 


2 
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where  n  8  and  s  is  the  sample  variance. 

Two  sources  ol'  error  in  the  measurements  of  p.  were  mentioned  in 

1 

Section  C.  One  is  introduced  by  the  motion  of  the  aircraft,  which  was 
shown  to  introduce  an  rm.s  error  into  the  doppler  measured  by  the  aircraft 
(ffp)  of  about  0.02  Hz.  The  rms  measured  doppler  averaged  over  all 
flights  was  approximately  0.1  Hz.  Therefore  the  expected  error  in  p 
is  0.2  for  this  source.  The  expected  error  in  p  is  then  0.2/  8,  or 
approximately  0.07. 

The  other  source  of  error  is  introduced  by  the  horizontal  bending 
of  ray  paths  by  irregularities,  which  was  shown  to  introduce  an  uncer¬ 
tainty  in  the  location  of  the  ray  path  of  about  7  km  near  the  reflection 
height  of  the  rays .  This  can  be  represented  as  an  error  on  the  ground 
ol  11  km  for  each  ray.  The  error  in  Ax  will  depend  on  how  the  devia¬ 
tions  of  the  two  rays  are  correlated.  At  very  small  (compared  to  the 
size  of  the  irregularities)  spacings  the  deviations  are  positively  cor¬ 
related,  decreasing  the  error  in  Ax.  At  very  large  spacings,  the  devia¬ 
tions  will  be  independent.  At  some  intermediate  spacings,  the  two  rays 
will  be  deflected  oppositely  by  the  two  sides  of  each  irregularity,  caus¬ 
ing  the  deviations  to  be  negatively  correlated,  thus  increasing  the  error 
in  Ax.  On  the  basis  that  the  deviations  are  uncorrelated,  the  expected 
error  in  p  is  represented  as  an  error  in  the  location  of  p(£x).  This 
error  is  approximately  7  km  (14  8).  The  effects  of  this  error  and 

the  error  due  to  aircraft  motion  are  shown  in  Fig.  7. 

The  data  of  Fig.  7  show  a  correlation  greater  than  0.8  for  the  dop¬ 
pler  measured  at  spacings  of  60  km  or  less.  At  spacings  greater  than  00 
km,  the  correlation  is  less  than  0.3.  This  indicates  that  the  irregular¬ 
ities  producing  the  major  component  of  the  observed  dopplers  are  of  limi¬ 
ted  extent,  and  that  the  experimental  configuration  covered  the  size  range 
of  the  major  irregular  structure.  At  a  receiver  separation  of  70  km,  the 
average  correlation  has  dropped  to  about  0.5.  This  gives  us  a  reasonable 
first  estimate  of  35  km  for  the  horizontal  extent  of  the  irregularities 
in  the  ionosphere.  This  scale  size  will  be  used  as  a  typical  value  in 
predicting  the  expected  effects  which  these  irregularities  are  expected 
to  produce  on  wave  fronts. 
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As  was  noted  above,  the  doppler  measured  by  the  aircraft  is  a  func¬ 


tion  not  only  of  the  time  rate  oa  change  of  ray  path  length  at  a  fixed 
point,  but  also  of  the  spatial  variation  of  refractive  index  caused  by 
the  irregularities.  The  effect  of  the  spatial  structure  of  the  irregu¬ 
larities  of  20  min  period,  identified  in  the  above  doppler  correlation 
study,  was  estimated  earlier  in  this  section  using  Sweeney's  angle-of- 
arrival  measurements,  and  was  found  to  introduce  only  a  small  unbiased 
error  term  in  the  estimate  of  correlation.  The  possibility  of  the  exis¬ 
tence  of  a  fixed  irregular  structure  (or  structure  varying  with  a  period 
much  greater  than  the  observed  20  min  period)  was  examined  by  correlating 
the  doppler  measured  on  two  pairs  of  consecutive  flights  over  the  path 
(spaced  in  time  by  about  30  min).  Doppler  correlations  were  computed 
from  the  data  obtained  on  each  pair  of  flights.  The  correlation  coeffi¬ 
cients  of  0.13  and  -0.05  for  the  two  pairs  of  flights  indicate  that  no 
fixed  or  semi-fixed  structure  was  present  with  a  magnitude  comparable  to 
that  of  the  irregularities  having  the  20  min  period  variations. 

A  more  accurate  measure  of  the  scale  size  of  the  irregularity  struc¬ 
ture  could  be  obtained  if  it  were  possible  to  make  prolonged  measurements 
at  fixed  receiver  separations.  If  multiple  spacings  were  employed,  some 
estimate  of  the  spatial  spectrum  of  irregularities  could  be  obtained. 

Such  a  program  could  prove  costly  because  of  the  large  spacings  needed 
and  the  high  stability  requirements  of  each  receiving  system.  For  the 
purposes  of  this  work,  the  approximate  scale  size  of  the  irregularity 
structure  determined  above  will  be  used  to  study  the  effects  to  be  ex¬ 
pected  in  the  presence  of  such  irregularities. 


E.  VERTICAL  SCALE  SIZE 

Having  determined  a  reasonable  estimate  of  the  strength  and  of  the 
effective  horizontal  scale  size  of  the  most  significant  irregularities 
for  the  oblique  path,  it  was  then  necessary  also  to  determine  the  approx¬ 
imate  vertical  extent  of  the  irregularities  in  order  to  model  their  ray- 
beinding  effects.  The  most  convenient  technique  for  probing  the  iono¬ 
sphere  at  various  heights,  for  a  path  joining  two  points,  is  to  make  use 
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oi  a  ran nu  01  radio  .frequencies  which  will  penetrate  to  various  heights 
into  the  ionosphere.  Figure  8  shows  the  geometry  of  ray  paths  of  various 
frequencies  between  two  points,  and  gives  the  values  of  ray  height  at 
different  points  along  the  path,  as  computed  by  raytracing  techniques 
(Croft  '69),  ior  an  ionosphere  typical  of  that  over  the  path  under  study. 
These  heights  are  used  in  estimating  the  vertical  scale  size  of  the  ir¬ 
regularities  . 

During  the  winter  and  spring  of  1969-1970  another  series  of  experi¬ 
ments  was  performed  over  the  same  path  used  in  the  previous  work,  with 
the  objective  of  determining  the  vertical  scale  size  of  the  irregulari¬ 
ties  previously  observed.  Only  the  Los  Banos  receiver  site  was  used  for 
these  measurements.  Transmissions  on  two  single-mode  frequencies  were 
made  from  Bearden  and  received  at  Los  Banos,  The  two  frequencies  were 
received  on  a  single  whip  antenna.  One  of  the  frequencies  was  also  re¬ 
ceived  on  a  whip  antenna  located  2.52  km  north  of  the  first  antenna,  for 
use  in  measuring  the  north-south  velocity  component  of  the  irregularities. 
The  outputs  of  the  three  receivers  were  recorded  using  the  same  phase¬ 
measuring  technique  used  previously  (see  Appendix  A).  During  the  period 
from  November  1969  to  April  1970,  measurements  were  made  one  or  two  days 
per  week.  Each  day,  from  two  to  six  hours  of  operation  were  possible 
(periods  when  single-mode  propagation  was  possible  and  no  irregularities 
were  observable  on  ionograms  over  the  path).  The  data  collected  in  the 
course  of  these  measurements  were  also  analyzed  to  determine  the  vertical 
and  north-south  velocity  components,  seasonal  occurrence,  and  signal-fad¬ 
ing  effects  oi  these  irregularities.  These  results  are  discussed  in  later 
chapters  of  this  report. 

Frequency  separations  in  the  range  of  100  kHz  to  6  MHz  were  used  to 
achieve  various  ray  spacings  in  the  ionosphere.  Doppler  correlations 
similar  to  those  employed  in  the  determination  of  horizontal  scale  size 
were  performed  on  pairs  of  frequencies  using  1.2  hr  segments  of  data.  A 
considerably  greater  quantity  of  data  were  available  for  this  work  (75 
segments  of  1.2  hr  each)  than  were  available  from  the  aircraft  measure¬ 
ments  . 

Similar  phase-path  variations  were  observed  at  all  frequencies  used 
in  these  measurements.  (Examples  of  the  phase-path  records  are  shown  in 


Chapter  III,  Figs.  17-20.)  Frequently,  however,  the  amplitude  of  the 
phase-path  changes  decreased  with  increasing  frequency.  This  information 
may  be  of  use  in  determining  the  mechanism  responsible  for  the  observed 
variations;  but  for  the  present,  only  the  vertical  scale  size  of  the  ir¬ 
regularities  is  of  interest.  The  form  of  the  correlation  computed  was 
slightly  modified  in  order  to  utilize  the  greater  quantity  of  data  avail¬ 
able  and  to  compensate  for  the  general  decrease  in  magnitude  with  frequency. 
The  denominator  used  in  calculating  the  correlation  coefficient  was  changed 
from  f j  as  in  Eq .  (7)  to  if?  .  f?  This  form  of  the  correlation  func- 
tion  has  the  property  that  two  functions  differing  only  by  a  constant  fac¬ 
tor  will  have  a  correlation  coefficient  of  1.0.  For  each  sequence  of  data, 
taken  at  a  fixed  frequency  spacing,  a  correlation  was  calculated  from  64 
doppler  measurements  during  the  1.2  hour  period.  In  Fig.  9,  the  correla¬ 
tions  are  plotted  against  the  frequency  separation  used.  On  the  assump¬ 
tion  (supported  in  the  next  chapter)  that  the  irregularities  have  their 
maximum  effect  near  the  maximum  altitude  of  the  rays,  spatial  separations 
are  also  assigned  to  the  frequency  axis  in  Fig.  9  by  using  an  average 
value  of  4  km  per  MHz.  The  error  in  this  average  spacing  as  determined 
by  ray tracings  in  several  representative  model  ionospheres  (see  Fig.  8) 
is  indicated  on  the  figure. 

Again,  a  smoothly  decreasing  correlation  is  observed  as  the  spacing 
increases.  The  average  vertical  distance  required  for  the  doppler  cor¬ 
relation  to  drop  to  0.5  is  about  15  km,  based  on  the  assumed  ray  separa¬ 
tions.  This  considerably  smaller  vertical  scale  size  of  the  irregulari¬ 
ties,  in  comparison  to  horizontal  scale  size,  will  be  of  importance  in 
determining  the  effects  produced  by  these  irregularities. 


F.  CONCLUSIONS  FROM  PHASE-PATH  MEASUREMENTS 

The  important  parameters  of  the  irregularities  observed  will  now  be 
summarized  in  a  form  which  will  allow  an  investigation  in  the  following 
chapters  of  the  expected  effect  of  these  irregularities  on  wave  front 
characteristics  and  antenna  performance.  All  of  the  observations  of  ir¬ 
regularities  during  the  course  of  this  work  were  made  under  undisturbed 
winter  daytime  conditions  over  a  2600  km  path  at  frequencies  for  which 
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Figure  9.  Doppler  correlation  coefficient  for  multiple-frequency 
phase-path  measurements. 
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only  the  one-hop  lower-rav  F-layer  mode  was  propagating.  No  data  were 
taken  during  periods  when  larger  irregularities  (which  are  visible  as 
kinks"  in  the  ionogram — Fig.  3b)  were  present.  The  conditions  represent 
the  most  promising  conditions  for  high-performance  HF  propagation,  and 
are  therefore  of  great  practical  interest.  Such  conditions  also  allow 
the  use  of  the  most  direct  measurement  methods.  Furthermore,  the  appli¬ 
cation  of  the  knowledge  gained  under  these  conditions  will  help  to  clar¬ 
ify  the  results  of  observations  made  under  more  complex  conditions. 

During  about  100  hours  of  observation  in  the  spring  of  1969,  and  in 
the  winter  and  spring  of  1969-1970,  phase-path  variations  having  a  dis¬ 
tinct  period  of  15  to  25  min  were  always  observed.  The  amplitude  of  these 
variations  ranged  from  about  8  A  to  75  X  peak  to  peak.  The  horizontal 
spatial  correlation  of  the  doppler  produced  by  these  variations  indicates 
a  correlation  of  0.5  for  receiver  spacings  of  about  70  km  on  the  ground. 
This  implies  a  horizontal  scale  size  for  the  irregularities  of  about  35 
km.  (If  scale  size  is  defined  for  a  correlation  of  1/e,  where  e  is 
the  base  for  natural  logarithms,  the  scale  size  would  be  about  43  km.) 

A  similar  calculation  of  vertical  scale  size  indicates  a  vertical  extent 
of  about  15  km,  assuming  the  effects  are  produced  at  the  maximum  height 
of  the  rays.  Also,  for  a  fixed  ionosphere,  a  decrease  in  magnitude  of 
the  phase-path  variation  with  increasing  radio  frequency  has  been  noted. 
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III.  MODEL  OF  THE  IONOSPHERIC  IRREGULARITIES 


The  p  lase-path-length  measurements  described  in  the  preceding 
chapter  were  used  to  determine  the  strength  and  spatial  extent  of  the 
irregularities  producing  the  measured  variations  in  phase-path  length. 
The  measurements  do  not,  however,  provide  a  complete  description  of  the 
received  wave  field.  In  fact,  the  measuring  technique  intentionally 
rejects  amplitude  fluctuations  and  small  phase  fluctuations  in  order  to 
provide  measurements  which  describe  the  causative  mechanism  (the  irregu¬ 
larities)  most  directly.  In  this  chapter,  irregularity  models  are  used 
to  determine  the  nature  of  the  total  wave  field  expected  in  the  presence 
of  the  weak,  medium  scale  irregularities. 

A  two-dimensional  representation  of  ray  bending,  as  shown  in  Fig. 
10,  will  be  used.  In  the  figure,  r  is  the  distance  from  the  scatter¬ 
ing  region  to  the  point  of  observation,  and  e  is  the  angular  deviation 
imposed  on  the  ray  by  the  irregularities.  This  representation  will  be 
applied  to  diffraction  in  both  the  vertical  and  horizontal  planes. 


A.  THE  RAY-THEORY  APPROXIMATION 

The  spatial  extent  of  an  irregularity  strongly  influences  its  effect 
on  waves  which  propagate  through  it .  In  particular,  the  relation  between 
the  transverse  scale  size  (Ap)  of  the  irregularities  (measured  trans¬ 
versely  to  the  ray  path)  and  the  length  of  the  first  Fresnel  zone  (F) 
determines  the  suitability  of  a  ray-theory  treatment  (AT  »  F)  or  a 
diffraction  theory  treatment  (F  »  Ap) •  If  it  is  assumed  that  parallel 
rays  are  incident  on  the  scattering  region,  then 


F  =  ,  (9) 

where  r  is  the  distance  from  the  receiver  to  the  scattering  region  and 
A  is  the  wavelength.  For  the  parameters  of  the  measurements  in  this 
work,  r  ~  1300  km  and  A  ~  12  m.  Therefore,  F  «  4  km.  Since  the 
scale  size  of  the  irregularities  in  question  was  determined  in  the  pre¬ 
vious  chapter  to  be  considerably  larger  than  4  km,  a  ray-theory  treat¬ 


ment  is  appropriate . 
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TILE  EFFECTS  OF  A  SINGLE  IRREGULARITY 


In  this  section  it  is  assumed  that  the  observed  variations  in  phase- 
path  length  are  caused  by  single  isolated  circular  irregularities.  The 
refractive  index  within  the  irregularity  is  given  by 


where  ^  is  the  ambient  (background)  refractive  index  of  the  ionosphere, 
_Npmax  is  the  refractive  index  change  at  the  center  of  the  irregularity 
(the  strength  of  the  irregularity)  and  z  is  distance  from  the  center 
oi'  the  irregularity. 

1 .  Nature  of  the  Wave  Field  Produced  by  a  Single  Irregularity 

Computed  ray  paths  through  a  single  circular  irregularity  are 
shown  in  Fig.  11  for  three  different  values  of  irregularity  strength. 
Figures  lla,  b  and  c  illustrate  the  effect  of  successively  stronger  ir¬ 
regularities.  The  scales  are  greatly  compressed  in  the  direction  of 
propagation  in  order  to  provide  a  convenient  format. 

The  angle  e  through  which  rays  are  bent  when  passing  through 
such  an  irregularity  is  a  function  of  the  phase  change,  cp,  produced  by 
the  irregularity  and  of  the  transverse  scale  size  A^,  of  the  irregular¬ 
ity.  The  wave  field  at  the  receiver  plane  is  determined  by  the  angular 
deviations  of  the  rays  and  by  the  distance  from  the  irregularity  to  the 
receiver.  If  the  phase  change  of  a  ray  passing  through  the  center  of  an 

irregularity  is  defined  as  cp  ,  the  nature  of  the  received  wave  field 
13  max 

is  determined  by  the  ratio  of  cp  to  A  / F.  If  cp  is  small  com- 

max  T  max 

pared  to  A^/F,  only  a  single  ray  will  arrive  at  a  given  point  on  the 

ground  (Fig.  lla).  However,  if  cp  is  large  compared  with  A^F, 

max  l 

rays  will  cross  before  reaching  the  receiver  (Fig.  11c),  thus  producing 
a  field  at  the  receiver  plane  which  may  be  composed  of  three  separate 
waves  having  slightly  different  arrival  angles. 

By  calculating  points  of  constant  phase-path  length  (i.e., 
fj.;ds=  constant)  along  each  ray  path,  the  shape  of  the  resulting  wave 
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fronts  may  be  determined.  By  computing  the  change  in  area  of  each  flux 

tube  (an  arbitrary  volume  bounded  by  rays),  the  resulting  distributions 

of  field  amplitude  may  be  determined  (with  some  restrictions  in  regions 

where  rays  are  focused).  Examples  of  the  results  of  such  calculations 

are  shown  in  Figs.  12  and  13  for  irregularities  which  produce  maximum 

phase  changes,  cp  ,  of  cp  ,  4rp  ,  and  6.5rp  ,  where  cp  is  chosen 
max  11  l  l 

arbitrarily  simply  to  illustrate  the  effects  of  varying  irregularity 
strengths.  Note  that  the  wave-front  plots,  Fig.  12,  could  be  plotted  on 
the  ray  tracings,  Fig.  11;  but  the  scale  of  interest  is  drastically  dif¬ 
ferent.  For  the  irregularities  we  are  studying,  the  total  ray-path 
length  is  over  2600  km,  while  the  maximum  wave-front  deviation  is  about 
1  km. 

In  the  presence  of  a  single  irregularity,  the  field  at  the  re¬ 
ceiver  will  be  composed  of  one  or  three  waves,  depending  on  the  strength 
and  position  of  the  irregularity.  If  multiple  waves  are  present,  the 
amplitudes  of  the  waves  will  be  similar,  differing  by  only  a  few  dB.  As 
the  strength  or  position  of  the  irregularity  changes  with  time,  deep 
fades  will  be  observed  as  the  relative  phase  of  the  similar  amplitude 
waves  changes  with  time. 

2 .  Relation  of  Variations  in  Phase-Path  Length  to  Variations  in 

Refractive  Index  and  Electron  Density 

The  two-dimensional  single-irregularity  model  just  described 
will  now  be  used  to  estimate  the  variations  in  refractive  index  and  elec¬ 
tron'  density  required  to  produce  the  experimentally  observed  variations 
in  phase-path  length.  The  calculations  ignore  the  ray  curvature  produced 
by  the  slowly  varying  background  refractive  index.  However,  it  is  shown 
in  Appendix  B  that  these  calculations  nevertheless  provide  a  valid  rela¬ 
tion  between  observed  changes  in  phase-path  length  and  changes  in  refrac¬ 
tive  index.  Since  the  simple  model  assumes  that  all  observed  phase-path 
variation  is  caused  by  one  individual  irregularity,  this  calculation  will 
set  an  upper  bound  on  the  electron-density  variation.  The  peak-to-peak 
variations  in  P  observed  experimentally  were  in  the  range  of  100  to 
1000  m.  Assuming  that  the  irregularities  are  roughly  circular  in  plan 
view  (as  was  assumed  in  the  raytracing),  the  length  of  the  ray  path  in 
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Figure  13.  Computed  received  amplitude  for  waves 

propagating  through  a  gaussian  irregularity. 
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the  irregularity  is  about  35  km.  Therefore,  the  ray  path  is  modified  by 
about  100  to  1000  in  in  a  distance  of  35  km.  From  the  defining  equation 
ot  the  ray-path  Length,  Eq .  (1),  it  is  seen  that  AP/P  is  proportional 
to  AiV u .  Therefore,  the  variations  of  refractive  index  are  in  the  range 
oi  0.3  per  cent  to  3.0  per  cent.  The  relation  between  the  electron  den¬ 
sity  variation,  yN/N,  and  the  refractive-index  variation,  Ap/p,  is  a 
function  of  the  background  refractive  index,  p. 

For  the  ionosphere,  neglecting  electron/ion  collisions  and  the 
effect  of  the  magnetic  field  (which  are  unimportant  for  this  calculation) 
the  refractive  index  is  given  by 


80.6  N 


where  N  =  electron  density  in  electrons  per  cubic  meter,  and  f  =  fre¬ 
quency  in  MHz. 

Figure  14  is  a  plot  of  the  per  unit  change  in  electron  density, 
AN/N,  required  to  produce  a  one  per  cent  change  in  refractive  index, 
plotted  against  the  background  refractive  index,  p.  For  the  path  and 
frequencies  under  investigation,  the  refractive  index  along  the  ray  path 
varies  from  1.0  to  a  minimum  of  about  0.7  at  the  maximum  height  of  the 
ray.  At  a  background  refractive  index  of  0.7,  a  one  per  cent  change  in 
refractive  index  requires  about  a  two  per  cent  change  in  electron  density. 

At  points  of  lower  height  along  the  ray  path,  the  refractive  index  is 
greater,  and  larger  changes  in  electron  density  are  required  to  produce 
the  same  modification  of  refractive  index  and  hence  of  the  ray  path  length. 
For  the  changes  in  refractive  index  observed  in  the  present  studies,  elec¬ 
tron  density  variations  of  0.6  per  cent  to  6.0  per  cent  are  required  at 
the  reflection  height.  Considerably  greater  changes  would  be  required 
at  other  points  on  the  ray  path. 

The  calculated  changes  in  electron  density  may  be  compared  with 
the  variations  observed  by  others.  A  great  deal  of  the  information  avail¬ 
able  on  ionospheric  irregularities  has  been  acquired  with  vertical-incidence 
pulse  sounders.  However,  ionospheric  dispersion  limits  the  time  delay 
(range)  resolution  of  measurements  made  with  these  sounders.  Also,  antennas 
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having  little  directivity  are  usually  employed  with  such  sounders,  caus¬ 
ing  their  sensitivity  to  decrease  for  spatially  small  irregularities  be¬ 
cause  of  the  spatial  averaging  effect  of  the  antenna.  According  to  Munro 
(’58),  the  lower  limit  for  detection  of  medium-seals  irregularities  with 
vertical-incidence  pulse  sounders  is  "a  few  per  cent"  of  the  background 
electron  density.  It  is  therefore  reasonable  to  assume  that  if  the  ir¬ 
regularities  most  prominent  in  our  measurements  are  located  near  the  ray 
reflection  height,  they  would  not  bo  seen  by  vertical-incidence  sounders, 
or  at  best  would  be  at  the  lower  limit  of  resolution  for  vertical-incidence 
techniques.  We  have  already  noted  the  difficulties  involved  in  observing 
such  irregularities  using  doppler  Fourier  analysis  methods,  and  we  have 
concluded  that  the  irregularities  we  observe  are  at  or  below  the  lower 
resolution  limit  of  this  method . 

The  use  of  satellite  signals  in  measuring  small  or  weak  irregu¬ 
larities  is  also  limited.  Satellite  measurements  indicate  variations  in 
the  total  electron  content  of  the  ionosphere  between  the  observation  point 
and  the  satellite,  and  hence  are  somewhat  less  sensitive  in  detecting  ir¬ 
regularities  of  limited  vertical  extent  than  are  ground-based  measurements. 
Also,  satellite  measurements  do  not  directly  yield  the  altitude  at  which 
irregularities  occur.  Titheridge  ('03)  observed  variations  in  total  elec¬ 
tron  content  in  the  range  of  0.05  per  cent  to  7.0  per  cent,  which  for  the 
size,  shape  and  altitude  (somewhat  below  the  F-layer  peak)  of  the  irregu¬ 
larities  we  observe,  would  require  variations  in  electron  density,  with 
reference  to  jthe  ambient  background,  of  0.5  per  cent  to  about  25  per  cent. 

It  is  likely  that  some  of  the  smaller  irregularities  observed  by  Titheridge 
were  the  same  type  observed  in  the  present  work.  However,  satellite  methods 
do  not  frequently  observe  irregularities  as  small  (horizontally)  as  those 
observed  here,  because  of  the  low  probability  that  a  ray  path  between  a 
satellite  and  a  ground  station  will  pass  through  a  given  irregularity. 

In  summary,  then,  it  appears  that  the  irregularities  whose  ex¬ 
istence  and  characteristics  are  proposed  from  the  two-dimensional,  single¬ 
irregularity  model  are  of  a  size  and  strength  range  not  readily  observed 
using  other  techniques.  Such  irregularities  are,  however,  very  significant 
in  their  effect  on  HF  oblique-incidence  radio  communication,  as  will  be 
seen  in  Chapter  IV. 


3 .  Modifications  nl'  the  Simple  Two-Dimensional  Model  to  Account 

lor  Background  Variation  in  p 

So  lar,  the  effects  of  the  curvature  of  the  ionospheric  ray 
paths  in  the  vertical  plane  (vertical  focusing),  arising  from  the  slowly 
varying  refractive  index  of  the  background  as  the  ray  travels  upward  have 
not  been  considered.  The  simple  model  of  Fig.  10  is  directly  applicable 
to  horizontal  (azimuthal)  ray  bending,  if  the  value  of  F  is  modified 
appropriately  for  rays  diverging  from  a  point  source.  In  this  case 
V  =  v/rX/2,  ignoring  a  minor  error  due  to  the  curvature  of  the  earth  and 
the  ionosphere.  However,  a  realistic  model  for  vertical  ray  beinding 
should  include  vertical  focusing  effects  caused  by  the  smooth  variation 
with  height  of  the  ionospheric  refractive  index. 

While  ray  bending  effects  may  be  determined  for  any  specific 
irregularity  distribution  using  computerized  raytracing  methods,  it  is 
of  greater  value  in  the  present  research  to  develop  an  analytical  descrip¬ 
tion  of  these  effects.  The  correction  for  ionospheric  focusing  in  the 
vertical  plane  may  be  made  by  determining  the  length  of  the  first  Fresnel 
zone  (F)  in  the  vertical  plane  for  the  region  where  the  ray  paths  are 
affected  by  irregularities. 

The  fundamental  definition  of  F  is  the  distance  along  a  wave 
front  at  which  two  ray  paths  to  a  receiving  point  differ  by  one-half 
wavelength  (Born  and  Wolf,  ’64).  The  parameter  F  therefore  describes 
the  effect  that  phase  changes  introduced  at  different  points  across  the 
wave  front  will  have  on  the  field  at  large  distances  from  the  phase¬ 
changing  region.  By  computing  wave  fronts  produced  by  radiation  from  a 
transmitter  on  the  ground,  using  rays  traced  through  a  typical  ionosphere, 
the  distance  F  may  be  determined  at  various  points  along  the  ray  paths. 
This  was  done  to  determine  the  value  of  F  near  the  reflection  height 
for  a  number  of  ranges  and  frequencies  typical  of  long-distance  HF  com¬ 
munication.  In  Fig.  15  the  ratio  of  the  actual  value  of  F  to  s/rA  at 
the  reflection  height  is  plotted  for  various  ranges  and  frequencies.  For 
the  conditions  of  our  measurements,  the  ratio  of  F  to  \/rA  is  less 
than  2.0.  The  use  of  the  actual  value  of  F  will  allow  the  inclusion, 
in  the  following  calculations,  of  the  focusing  effects  produced  by  the 
smooth  variation  with  height  of  refractive  index  in  the  ionosphere. 
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Figure  15.  Computed  Fresnel  zone  length  at  the  reflection  height  for 
rays  in  the  ionosphere. 
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C.  THE  EFFECTS  OF  MULTIPLE  IRREGULARITIES  IN  THREE  DIMENSIONS 

The  two-dimensional  representation  of  propagation  and  as  affected 
by  irregularities  (as  in  Fig.  10)  will  now  be  applied  to  both  the  hori¬ 
zontal  and  vertical  axes  of  the  irregularities,  including  the  corrections 
for  focusing  determined  in  the  previous  section.  The  effects  of  multiple 
irregularities  will  also  be  included. 

1 .  The  Effect  of  the  Presence  of  Multiple  Irregularities  on  the 

Calculated  Change  in  Electron  Density 

The  ray tracings  shown  in  Fig.  11  are  the  results  of  raytracing 
through  a  single  irregularity.  There  may  be  more  than  a  single  irregu¬ 
larity  that  influences  wave  field  at  the  receiver.  However,  the  number 
of  such  irregularities  is  probably  not  large,  as  the  region  where  the 
ray  path  is  most  susceptible  to  irregularities  of  a  given  percentage  de¬ 
viation  in  electron  density  is  limited  to  the  region  near  the  maximum 
height  of  the  ray.  In  the  presence  of  an  isolated  two-dimensional  gaus- 
sian  irregularity,  the  number  of  rays  at  the  receiver  will  be  one  or 
three,  depending  on  the  strength  and  size  of  the  irregularity  (see  Fig. 
11).  In  the  presence  of  more  than  one  irregularity,  a  larger  number  of 
rays  is  expected  to  reach  the  receiver. 

Chernov  ('60)  calculates  the  mean  squared  phase  fluctuation 

2 

cp  ,  due  to  propagation  through  an  irregular  medium  described  by  a  spa¬ 
tially  stationary  gaussian  correlation  of  refractive  index  with  distance, 
as 

q)2  =  n/tT  (Ap)2  k2ATL  ,  (11) 

Li 

where 

Ap  =  the  change  in  refractive  index, 
k  =  2  jt/A, 

L  =  the  length  of  the  scattering  region, 

At  =  the  longitudinal  scale  size  of  the  irregularity. 

Li 
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From  the  measurements  made  in  this  work,  cp  may  be  calculated 
Using  this  estimate  in  (11),  (Al-O^  can  be  estimated,  assuming  that  the 
longitudinal  and  transverse  scale  sizes  of  the  irregularities  are  known. 
For  the  single-irregularity  model  used  previously,  it  was  assumed  that 
A,p  =  =  L.  While  these  assumptions  cannot  bo  verified  with  the  experi¬ 

mental  geometry  used,  they  do  however  produce  reasonable  values  for  Ap. 
If  more  than  one  irregularity  affects  a  single  ray  path  (L  >AT),  or 
if  the  irregularities  are  elongated  in  the  direction  of  propagation 
(A^  >  Aj.),  smaller  values  of  Ap  result.  This  would  be  acceptable  in 
terms  of  known  ionospheric  characteristics. 

2 .  The  Range  of  Phase  Changes  Produced  by  the  Observed 

Irregularities 

Ray  bending  effects  caused  by  irregularities  are  determined 
mainly  by  the  phase  change  produced  by  the  irregularities  and  the  trans¬ 
verse  (to  the  ray  path)  scale  size  A^  of  the  irregularities  (horizontal 
and  vertical).  A  variety  of  data  taken  during  the  experiments  of  this 
work,  which  is  representative  of  the  range  of  phase-path  variations  ob¬ 
served,  is  presented  here.  Figure  16  is  an  example  of  data  representing 
quiet  conditions.  Phase-path  length  changes  are  expressed  in  wavelengths 
Tne  record  shows  a  maximum  change  in  phase-path  length  of  about  20  wave¬ 
lengths  peak-to-peak .  The  trace  is  smooth,  indicating  that  a  single  ray 
was  dominating  the  measurement  throughout  the  period  of  the  observation. 
Figure  17  is  a  similar  record,  but  with  a  greater  observed  deviation  in 
phase-path  length.  The  record  is  still  dominated  by  a  single  ray  through 
out . 

Some  other  effects  occasionally  appear  in  the  records.  Figure 
18  shows  the  effect  of  alternate  capture  of  the  measurement  by  two  rays 
having  different  doppler  shifts.  The  record  which  would  be  produced  by 
either  of  these  rays  separately  is  shown  by  the  dotted  lines.  In  this 
record,  the  measurement  is  still  yielding  a  useful  description  of  the 
changes  in  phase-path  length  for  the  rays.  Finally,  Fig.  19  is  an  exam¬ 
ple  of  a  record  taken  under  conditions  in  which  the  usefulness  of  the 
measuring  technique  has  broken  down  because  of  the  presence  of  a  number 
of  rays  of  different  doppler  shifts,  each  of  which  dominates  the  measure¬ 
ment  for  short  periods  of  time. 
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PHASE- PATH  LENGTH,  P  (wavelengths) 
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Figure  16. 


Phase-path  length  measurement  for  very  quiet 
ionospheric  conditions. 
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Figure  17. 


Phase-path  length  measurement  for  fairly  quiet 
ionospheric  conditions. 
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Figure  18.  Phase -path  length  measurement 
alternately  captured  by  two 
rays  of  comparable  amplitude. 
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Figure  19.  Phase-path  length  measurement 
for  disturbed  ionospheric 
conditions . 
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3 .  The  Effects  of  the  Observed  Phase  Changes 

With  the  observed  range  of  phase  variations,  cp2,  in  mind, 
the  theoretically  predicted  wave-front  effects  in  the  presence  of  multi¬ 
ple  irregularities  can  be  determined.  Chernov  calculates  the  mean  squared 
angular  ray  deviations  as  a  function  of  cp2  and  of  A—  >  the  trans¬ 

verse  scale  size  of  the  J  regularities . 


2  1  /  A  \  2 

e  =  ~2  [— h5 

rt  V'T/ 


Rays  emerging  from  the  scattering  region  are  unchanged  in  am¬ 
plitude  or  position,  but  have  a  change  in  angle  imposed  upon  them.  As 
the  wave  propagates  through  space  to  the  receiving  location,  the  rays 
converge  or  diverge  (Fig.  11).  Also,  as  the  wave  propagates,  the  ampli¬ 
tude  at  various  positions  across  the  wave  will  change  due  to  the  conver¬ 
gence  and  divergence  of  the  rays.  Manning  ('68)  describes  two  classes 
of  wave  fields  which  result  from  this  process.  The  distinction  is  based 
upon  the  number  of  rays  which  are  expected  to  arrive  at  one  point  on  the 
ground.  For  weak,  large  or  nearby  irregularities,  only  a  single  ray  will 
be  expected  to  arrive  at  any  one  point  on  the  ground.  For  stronger,  smal¬ 
ler  or  more  distant  irregularities,  rays  will  cross  before  reaching  the 
plane  of  the  receiver,  and  multiple  rays  will  arrive  at  least  some  points 
on  the  ground.  Manning  calculates  the  expected  number  of  rays  at  a  given 
point  (for  the  two-dimen  ional  geometry)  as 


r  i  sfe  F  2 

EB3  =  i  + 


3/2  r~ 

where  E[*]  denotes  expectation.  The  factor  (rt  /V2)  =»  4.  So,  if 
cp2  «  16  (Ap/F  /»  it  is  probable  that  only  one  ray  will  arrive  at  a  given 
point  on  the  ground;  whereas  if  cp2  »  IBCA^/f/,  more  than  one  ray  can 
arrive  at  the  same  pointT 
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For  the  case  where  the  scattering  region  contains  irregulari¬ 
ties  which  have  variations  in  two  dimensions  transverse  to  the  ray  path, 
Manning  further  shows  that 

E[R]  =  E[Ry]  •  E[Rh]  ,  (14) 

where  E[Ry]  =  the  expected  number  of  rays  produced  by  vertical  bending, 

and  E  [R  1  =  the  expected  number  of  rays  produced  by  horizontal  bending. 

As  was  previously  pointed  out,  the  effective  value  of  F  is 

modified  in  the  case  of  vertical  ray  deviations  by  the  focusing  effect 

of  the  smooth  variation  in  |i  with  height  in  the  ionosphere.  By  using 

the  corrected  value  of  F  for  vertical  bending  as  plotted  in  Fig.  15, 

and  F  =  '\/rA/2  for  horizontal  bending,  the  calculations  of  the  expected 

number  of  rays,  as  well  as  the  signal  strength  calculations  (Fig.  13), 

can  be  corrected  for  the  effect  of  vertical  bending  in  the  ionosphere. 

For  the  observations  made  in  this  work,  cp2  was  found  to  vary 

over  the  range  of  approximately  100  to  10,000  radians  squared.  For  a 

horizontal  scale  size  of  45  km  and  a  vertical  scale  size  of  20  km,  E[R], 

E [R  1  and  E [R  ]  are  plotted  against  in  Fig.  20.  It  can  be  seen 

V  H 

that  the  major  contribution  of  the  ionosphere  to  the  multiple  rays  arriv¬ 
ing  at  the  receiver  is  refraction  in  the  vertical  plane,  and  that  for  al? 
but  the  quietest  conditions  seen  during  our  observations,  E [R]  is 
greater  than  1,  but  usually  less  than  10.  The  value  of  E[r]  from  Fig. 

20  can  be  regarded  as  an  upper  limit  on  the  expected  number  of  rays  reach¬ 
ing  the  receiver.  Therefore  it  is  expected  that  wave  fields  such  as  those 
of  Figs,  lib  and  11c  will  be  common,  and  in  the  presence  of  more  than  a 
single  irregularity,  somewhat  more  complex. 

The  total  amplitude  and  phase  of  the  wave  field  at  a  point  on 
the  ground  can  now  be  described  in  terms  of  the  summation  of  R  waves. 

The  relative  amplitudes  of  the  waves  will  be  in  the  range  shown  in  Fig. 

13.  The  resulting  total  amplitude  and  phase  fluctuations  are  described 
in  Chapter  IV. 


55  SEL-70-037 


EXPECTED  NUMBER  OF  RAYS, 


D. 


EFFECT  OF  THE  EARTH’S  MAGNETIC  FIELD 


Thus  far  the  effect  of  the  earth's  magnetic  field  on  oblique-inci¬ 
dence  HF  ionospheric  propagation  has  been  ignored.  This  effect  may  be 
described  either  in  terms  of  the  rotation  c-i  the  plane  of  polarization 
of  a  linearly  polarized  wave  (Faraday  rotation)  or  in  terms  of  splitting 
of  the  wave  incident  on  the  ionosphere  into  two  components,  each  of  which 
exhibits  one  of  two  possible  characteristic  polarizations  for  a  wave  of 
that  direction  and  frequency  in  the  ionosphere,  and  propagates  over  a 
slightly  different  path.  The  field  produced  by  the  combination  of  these 
two  component  waves  at  any  point  on  the  path  produces  the  linearly  polar¬ 
ized  wave  described  by  the  first  method. 

The  concept  of  Faraday  rotation  is  frequently  valid  for  oblique- 
incidence  ionospheric  propagation.  As  a  wave  propagates  through  the 
ionosphere  its  plane  of  polarization  (assuming  that  the  transmitted  wave 
is  plane-polarized)  rotates  slowly  with  distance  in  the  refracting  region 
of  the  ionosphere.  This  rotation  causes  the  wave  arriving  at  the  receiver 
to  have  a  plane  of  polarization  determined  by  the  amount  of  Faraday  rota¬ 
tion  which  occurred  in  the  ionosphere.  If  the  receiving  antenna  is  lin¬ 
early  polarized,  the  received  signal  strength  will  vary  as  the  plane  of 
polarization  of  the  received  wave  rotates.  For  the  path  used  in  our  ex¬ 
periments,  the  total  number  of  rotations  of  the  plane  of  polarization 
along  the  path  of  the  ray  is  calculated  by  raytracing  (Croft  ’69)  to  be 
about  20  turns  under  typical  daytime  conditions. 

It  is  not  possible  to  develop  a  simple  relation  between  variations 
in  the  total  phase-path  length  and  changes  in  the  number  of  Faraday  turns, 
for  a  ray  in  the  ionosphere.  Therefore,  various  assumed  ionospheric  elec¬ 
tron  density  profiles,  representative  of  various  times  of  day,  were  used 
for  tracing  rays,  while  calculating  both  the  total  phase  path  and  the  total 
number  of  Faraday  turns.  A  variation  of  50  to  100  wavelengths  in  P  pro¬ 
duced  a  change  of  one  turn  in  Faraday  rotation.  This  implies  that  for 
smooth  changes  in  the  ionosphere,  only  one  or  two  nulls  in  received  signal 
strength  should  occur  for  each  50  to  100  wavelengths  of  observed  change  in 
P. 

When  considering  the  effects  of  irregularities,  it  is  more  revealing 
to  consider  the  model  consisting  of  two  component  rays  propagating  over 


slightly  different  paths,  than  it  is  to  use  the  concept  of  Faraday 
rotation.  The  two  rays  are  customarily  denoted  as  ordinary  (0)  and 
extraordinary  (X).  If  these  rays  have  phase-path  lengths  Pq  and  P  , 
a  change  in  the  difference  between  the  two  phase-path  lengths  (P  -  P  ) 
of  one  wavelength  will  cause  one  Faraday-turn  null  at  the  receiver.  By 
tracing  rays  using  the  ESSA  three-dimensional  ray-tracing  program  (Jones 
'681,  it  was  determined  that  for  the  conditions  of  our  measurements,  the 
0  and  X  components  will  be  separated  vertically  by  a  maximum  of  0.5  km 
to  1.0  km  at  their  greatest  height  in  the  ionosphere.  Laterally,  they 
will  be  separated  by  no  more  than  200  m.  These  numbers  represent  the  ray 
spacings  under  undisturbed  ionospheric  conditions,  i.e.,  with  no  irregu¬ 
larities  present.  If  an  irregularity  causes  the  phase-path  length  of  one 
of  the  rays  to  change  with  respect  to  that  of  the  other  ray  by  one  wave¬ 
length,  one  null  will  be  observed  at  the  receiver.  Based  on  the  measured 
irregularity  scale  sizes,  and  assuming  gaussian  irregularities,  the  changes 
in  phase-path  length  for  two  rays  spaced  by  about  1  km  vertically  are 
highly  correlated — of  the  order  of  0.98.  Therefore,  for  irregularities 
such  as  those  studied  and  modeled  in  this  work,  about  one  Faraday  null  is 
expected  for  each  50  wavelengths  change  in  P  due  to  irregularities. 

Based  on  the  above  calculations,  then,  both  irregularities  and  smooth 
variations  in  the  ionosphere  are  expected  to  produce  variations  in  the 
received  polarization  which  will  cause  one  receiver  null  for  each  50  to 
100  wavelengths  change  in  P.  This  sets  an  upper  limit  to  the  number  of 
such  nulls  which  can  be  explained  on  the  basis  of  the  measured  irregular¬ 
ity  sizes. 
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IV.  SIGNAL- STRENGTH  FLUCTUATIONS 


In  Chapter  II  a  series  of  measurements  of  the  phase  characteristics 
of  HF  obliquely  propagated  ionospheric  waves  was  described.  These  mea¬ 
surements  allowed  a  determination  of  the  size  and  strength  range  of  ir¬ 
regularities  which  are  present  in  the  ionosphere  under  quiet  conditions. 
In  Chapter  III  the  effects  which  would  be  anticipated  for  waves  propagat¬ 
ing  through  such  irregularities  were  examined.  The  expected  number  of 
rays  which  would  be  produced  by  these  irregularities  was  determined  and 
the  expected  relative  magnitudes  were  found.  In  the  present  chapter, 
the  amplitude  variations  which  are  observed  under  quiet  ionospheric  con¬ 
ditions  are  examined  to  verify  that  observed  amplitude  fluctuations  are 
explained  in  terms  of  interference  of  rays.  In  addition,  other  possible 
causes  of  the  observed  amplitude  fluctuations  are  discussed. 

The  envelope-amplitude  fluctuations  of  a  received  signal  can  be  de¬ 
scribed  by  their  statistical  distribution.  The  probability  density  func¬ 
tion  for  the  envelope-amplitude,  $(E),  will  be  calculated  for  various 
assumptions  about  the  received  waves.  These  calculated  distributions 
will  be  compared  with  histograms  of  observed  envelope -amplitude  fluctua¬ 
tions.  The  histograms  are  plotted  so  that  the  densities  of  the  measured 
values  of  envelope-amplitude,  W(E),  are  directly  comparable  with  the 
computed  probability-density  functions.  The  mean  envelope  amplitude  E 
is  used  as  a  scale  factor  in  plotting  both  the  computed  and  observed  en¬ 
velope-amplitude  distributions. 

In  addition  to  the  density  functions  themselves,  the  normalized 

2 

variance  of  the  envelope  amplitude,  a  ,  will  be  calculated  for  both 

2 

predicted  and  observed  amplitude  fluctuations.  The  value  of  a  is 
given  by 


A  second  statistic,  cp^,  is  also  calculated  for  each  of  the  assumed 
forms  of  the  received  wave  field.  Here  cp  is  the  difference  between 
the  phase  of  the  received  signal  at  any  instant  of  time  and  the  phase 
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of  the  strongest  of  the  signal's  component  waves.  The  values  of  o  and 
may  be  used  to  determine  the  efficiency  of  receiving  antennas  under 
certain  conditions.  This  subject  will  be  discussed  briefly  in  the  next 
chapter . 


A.  MULTIPLE  WEAK-SCATTERING  INTERPRETATION  OF  AMPLITUDE  FLUCTUATIONS 

Booker,  Ratcliffe  and  Shinn  ('50)  applied  diffraction  theory  to  the 
ionospheric  propagation  problem  by  describing  scattering  in  the  ionosphere 
in  terms  of  a  scattering  region  composed  of  a  thin  screen  of  small,  weak 
phase-changing  irregularities.  The  signal  which  emerges  from  such  a  scat¬ 
tering  region  can  be  characterized  in  terms  of  a  specular  component  and 
a  cone  of  scattered  rays.  The  "coherence  ratio",  p,  is  used  to  describe 
the  ratio  of  the  energy  of  the  specular  component  to  the  energy  of  the 
scattered  rays.  The  width  of  the  cone  of  scattered  rays  is  determined 
by  the  size  of  the  irregularities  in  relation  to  the  wavelength  employed. 
The  cone  of  rays  is  made  up  of  a  large  number  of  rays  having  random  am¬ 
plitudes  and  phases.  If  the  number  of  scattered  rays  reaching  the  re¬ 
ceiver  is  statistically  large,  the  expected  envelope -amplitude  distribu¬ 
tions  for  the  sum  of  a  sine  wave  and  narrow-band  random  noise,  as  derived 
by  Rice  ('45),  may  be  used  to  predict  analytically  the  amplitude  fluctua¬ 
tions  produced  by  the  specular  and  scattered  components  of  ionospherically 
propagated  waves.  Rice  derives  the  probability  density  for  the  envelope- 
amplitude  of  such  a  signal  as 


~2  eXP 


2  2 
E  +  EZ, 
a 


where  (in  the  notation  used  in  the  present  report): 

E  =  the  envelope  amplitude  of  the  received  signal, 

E  =  the  amplitude  of  the  specular  component, 

S 

E  =  the  amplitude  of  the  scattered  component, 
n 


and 


I  is  the  zero-order  Bessel  function  of  imaginary  argument. 

Curves  of  0(E),  using  the  normalizations  employed  by  Rice,  are  shown 
for  different  values  of  the  parameter  s  =  in  Fig.  21.  Before  an 

experimental  amplitude  density  distribution  may  be  compared  with  these 
curves,  the  parameter  (3  must  be  determined.  (3  is  related  to  the  mea- 
sureable  quantity  E2/(E)2,  the  ratio  of  the  mean  squared  amplitude  to 
the  squared  mean  of  the  amplitude,  by  the  relation 


where  I  and  I„  are  Bessel  functions  of  imaginary  argument.  This 

°  1  -o  -  2 

function  is  plotted  in  Fig.  22.  EV(E)  must  lie  in  the  range 

1  <  E2/(E)2  <  4/ fl. 

Experimentally  determined  amplitude  distributions  may  be  compared 
with  the  Rice  distributions  in  order  to  determine  the  applicability  of 
the  multiple  weak-scattering  theory  proposed  by  Booker,  Ratcliffe  and 
Shinn.  Because  of  the  necessity  for  calculating  (3  before  plotting  the 
experimental  distributions,  Rice's  method  of  displaying  amplitude  dis¬ 
tributions  (Fig.  21)  is  not  suitable  for  conditions  where  Eq .  (15)  may 
not  apply.  Therefore,  the  Rice  distributions  arc  replotted  in  Fig.  23 
for  various  values  of  (3,  using  the  mean  amplitude  E  as  a  scale  fac¬ 
tor.  These  graphs  were  computed  using  a  random-number  generator  and  the 
Monte  Carlo  method.  This  computation  method  will  be  used  later  in  con¬ 
nection  with  other  models  of  amplitude  fading.  Note  that  for  (3  equal 
to  0  (a  random  component  only),  the  distribution  is  the  Rayleigh  distri¬ 
bution  . 

The  results  of  applying  the  multiple  weak-scattering  theory  to  ob¬ 
lique  F-layer  ionospheric  propagation  have  not  provided  a  consistent 
picture  of  the  ionospheric  mechanism  responsible  for  the  observed  signal 
fluctuations.  Al'pert  ('63)  discusses  some  of  the  results  of  such  mea¬ 
surements  and  some  of  the  problems  encountered  in  their  interpretation. 
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Figure  23. 


The  Rice  amplitude  distribution  computed  using  the  Monte 
Carlo  method  and  plotted  with  the  mean  envelope  amplitude 
IT  as  a  scale  factor. 

(a)  3  =  0,  a2  =  0.272. 

(b)  0  =  2,  a2  =  0.101. 

(c)  3  =  3,  a2  =  0.054. 

(d)  0  =  5,  a2  =  0.024. 
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Measurements  ol  (3  are  reported  in  the  range  from  0  to  10,  with  usual 
values  in  the  range  of  2  to  4.  No  numerical  evaluation  of  how  well  the 
observed  data  fit  the  predicted  results  of  the  multiple  weak-scattering 
theory  is  presented.  However,  Al'pert  notes  that  in  about  30  per  cent 
of  the  observations,  the  value  of  E^/(E)  exceeds  4/n,  the  maximum 
value  permitted  by  the  theory.  Also,  in  about  20  per  cent  of  the  obser¬ 
vations,  distributions  with  more  than  one  maximum  are  observed.  These 
discrepancies  lead  Al’pert  to  conclude  that  some  other  mechanism  must  be 
contributing  to  the  observed  signal-strength  fluctuations. 

The  same  ionospheric  scattering  model  used  in  the  amplitude  fluctua¬ 
tion  studies  discussed  above  has  been  applied  by  others  in  analyzing  the 
phase  fluctuations  observed  on  pairs  of  spaced  antennas.  Whale  and  Gar¬ 
diner  (’66)  derive  unique  probability  distributions  for  the  difference 
in  phase  measured  at  two  spaced  antennas  as  a  function  of  (3  and  the 
spatial  autocorrelation  function  of  the  scattered  component  alone.  Whale 
and  Boys  (’68)  and  Boys  (’68)  use  this  method  and  find  for  one-hop  propa¬ 
gation,  that  p  is  usually  quite  high  with  measured  values  of  10  or  more, 
indicating  that  the  received  wave  is  highly  specular.  The  results  of 
these  phase  measurements  are  not  compatible  with  those  commonly  obtained 
from  amplitude  measurements  using  the  Rice  theory. 

B.  INTERPRETATION  ON  THE  BASIS  OF  LARGE  IRREGULARITIES 

The  signal-strength  fluctuations  which  are  expected  in  the  presence 
of  large  irregularities  will  now  be  examined  in  terms  of  the  amplitude 
distributions  which  could  be  produced  by  such  irregularities.  Computed 
envelope-amplitude  probability-density  distributions  3>(E/E)  will  be 
used  to  gain  some  insights  into  the  effect  of  varying  parameters  on  the 
shapes  of  these  distributions.  Also,  the  normalized  variance  of  the  am¬ 
plitude  fluctuations  will  be  computed  for  comparison  with  the  known  lim¬ 
its  for  this  statistic  under  the  random-noise  (multiple  weak-scattering) 
representation.  Measurements  of  the  density  distribution  of  observed 
amplitude  fluctuations,  W(E/E),  will  be  compared  with  the  computed  pro¬ 
bability-density  distributions.  Four  different  possible  sources  of  sig¬ 
nal  amplitude  variation  will  be  discussed.  These  are  focusing,  multipath 
interference,  polarization  mismatch,  and  absorption. 


1 .  Focusing 


The  spatial  variations  in  refractive  index  produced  by  iono¬ 
spheric  irregularities  cause  changes  in  the  convergence  or  divergence  of 
rays.  This  in  turn  causes  variations  in  the  energy  density  between  adja¬ 
cent  rays.  This  change  in  energy  density  will  produce  variations  in  the 
observed  strength  of  a  given  wave  at  a  distance  from  the  refracting  region. 
As  can  be  seen  from  the  raytracings  in  Fig.  11,  the  variations  in  amplitude 
increase  with  distance  from  the  scattering  region.  These  variations  are 
in  the  range  of  a  few  decibels  for  the  irregularities  which  we  have  pro¬ 
posed.  (See  Fig.  13.)  In  cases  where  only  a  single  ray  is  present  at 
any  point  in  the  receiver  plane,  we  would  therefore  expect  to  see  signal- 
strength  variations  of  this  order  of  magnitude.  In  the  case  of  multiple 
rays,  it  is  noteworthy  that  the  difference  in  amplitudes  for  the  various 
components  is  again  only  a  few  d&.  This  fact  will  be  of  considerable  im¬ 
portance  in  considering  multipath  interference. 

If  changes  in  received  amplitude  are  attributable  only  to  varia¬ 
tions  in  the  strength  of  a  single  component  due  to  focusing,  the  expected 
fluctuations  in  signal  strength  would  be  a  few  dB.  This  is  observed  at 
times  on  HF  paths.  It  has  also  been  noted  in  ground  backscatter  radar 
observations.  Figure  24a  is  a  record  of  signal-strength  variations  of 
the  type  observed  over  the  path  used  for  the  phase  measurements  of  this 
report.  However,  signal-strength  fluctuations  as  small  as  these  are  much 
more  the  exception  than  the  rule.  Figures  24b  and  24c  show  examples  of 
more  typical  signal-strength  records  over  this  path.  Variations  such  as 
these  cannot  be  explained  on  the  basis  of  focusing  due  to  irregularities. 

2 .  Multipath  Interference 

Inasmuch  as  the  irregularity  model  developed  in  Chapter  III 
predicts  that  because  of  the  effect  of  even  a  single  irregularity,  the 
expected  number  of  rays  at  the  receiver  location  will  frequently  be 
greater  than  one,  the  effects  on  total  signal  strength  which  will  be 
produced  by  the  summation  of  several  rays  must  be  considered.  As  can 
be  seen  from  the  amplitude  plots  in  Fig.  13,  the  amplitudes  of  multiple 
rays  are  not  very  different  from  one  another.  Considering  the  vector 
addition  of  rays,  as  was  done  in  Chapter  II,  it  can  be  seen  that  the 
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Figure  24.  Observed  signal-strength  fluctuations  for  the  Arkansas-- 
Los  Banos  path. 

(a)  Little  fading. 

(b)  Moderate  fading. 

(c)  Severe  fading. 
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possible  range  of  amplitudes  for  the  sum  will  be  bounded  by  the  summation 
of  the  ray  amplitudes  as  an  upper  limit,  and  by  the  difference  between 
the  largest  ray  and  the  summation  of  the  remaining  rays  as  a  lower  limit. 

If  the  lower  limit  is  negative,  then  all  values  from  zero  to  the  upper 
limit  are  possible  (as  we  are  interested  in  the  envelope  amplitude).  The 
phases  of  the  multiple  rays  vary  sufficiently  that  in  performing  the  vec¬ 
tor  summations,  it  will  be  assumed  that  the  phase  of  each  ray  is  uniformly 
distributed  in  the  range  of  zero  to  2it.  This  assumption  is  usually  valid 
for  the  measurement  periods  used  in  this  work,  which  span  a  number  of  sig¬ 
nal  minima.  Computed  amplitude  probability-density  distributions  calcula¬ 
ted  for  two  rays  and  three  rays  having  amplitude  differences  of  a  few  dB 
(such  as  we  expect  to  encounter  based  on  the  work  of  the  previous  chapter) 
are  shown  in  Figs.  25a  and  b.  Note  that  these  amplitude  distributions  are 
quite  different  in  appearance  from  those  calculated  using  the  random-noise 
representation  for  a  large  number  of  small  scatterers  and  a  strong  specu¬ 
lar  component  (see  Fig.  23). 

3 .  Polarization  Mismatch 

As  was  pointed  out  in  Chapter  III,  the  ionosphere  is  birsfrin- 
gent  because  of  the  presence  of  the  earth's  magnetic  field.  Therefore, 
if  a  linearly  polarized  wave  is  incident  on  the  ionosphere,  it  will  be 
split  into  two  components  of  differing  elliptical  polarization,  and  each 
component  will  travel  a  slightly  different  path  through  the  ionosphere. 

Upon  emerging  from  the  ionosphere,  these  two  components  (the  ordinary 
and  extraordinary  rays)  will  add  to  produce  some  resultant  polarization 
which  is  a  function  of  the  difference  in  path  lengths  for  the  components, 
and  of  the  difference  in  strength  of  the  components.  At  the  end  of  Chap¬ 
ter  III,  it  was  shown  that  very  little  differential  in  path  length  or 
signal  strength  between  the  ordinary  and  extraordinary  rays  was  expected 
to  be  introduced  by  ionospheric  irregularities.  If  this  is  the  case, 
the  emerging  waves  will  sum  to  form  a  linearly  polarized  wave  (still 
assuming  that  the  incident  wave  was  linearly  polarized).  However,  the 
plane  of  polarization  will  no  longer  be  the  same  as  that  of  the  incident 
wave;  and  furthermore,  it  will  be  time-varying.  If  a  circularly  polar¬ 
ized  antenna  is  used  to  receive  such  a  wave,  its  response  will  be  constant, 
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Computed  amplitude  distribution  for  the  random  vector 
sum  of  rays. 

(a)  Two  rays,  amplitudes  =  1.0,  0.7. 
a2  =  0.175,  cp7  =  0.53. 

(b)  Three  rays,  amplitudes  =  1.0,  0.7,  0.5.  - 

a2  =  0.200,  cp7  =  0.80. 
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regardless  of  the  plane  of  polarization  of  the  wave.  In  practice,  it  is 
difficult  to  construct  and  test  an  antenna  designed  to  be  sensitive  to 
circular  polarization  for  the  frequencies  and  low  radiation  angles  used 
for  long-range  HF  communication.  If  a  linearly  polarized  receiving  an¬ 
tenna  is  used,  a  polarization  mismatch  will  result,  producing  variations 
in  response  depending  upon  the  instantaneous  plane  of  polarization  of  the 
received  wave.  The  signal  amplitude  from  such  an  antenna  will  be  propor¬ 
tional  to  cos  a,  where  a  is  the  angle  between  the  planes  of  polariza¬ 
tion  of  the  incident  wave  and  the  antenna.  The  two  components  (0  and 
X)  arriving  at  a  linearly  polarized  receiving  antenna  may  be  represented 
by  two  linearly  polarized  waves  having  amplitudes  equal  to  the  fraction 
of  each  wave's  energy  which  lies  in  the  plane  of  polarization  of  the  re¬ 
ceiving  antenna.  The  amplitudes  will  remain  fixed  over  long  periods  of 
time,  since  the  polarization  of  the  received  components  is  determined 
primarily  by  the  ray  direction  and  the  earth’s  magnetic  field.  Figure  26 
shows  the  calculated  amplitude  distributions  which  would  be  expected  for 
two  waves  of  differing  relative  amplitudes.  The  distributions  show  sharp 
upper  and  lower  limits  in  amplitude,  corresponding  to  the  sum  and  differ¬ 
ence  of  the  two  rays . 

4.  Absorption 

Waves  propagating  in  the  ionosphere  also  suffer  losses  due  to 
absorption  of  energy  associated  with  particle  collisions.  With  the  excep¬ 
tion  of  periods  of  severe  solar  storm  activity,  the  loss  due  to  absorption 
is  very  small  for  frequencies  above  a  few  megacycles.  Variations  in  loss 
from  this  mechanism  are  neglected  in  this  work. 

C.  SURVEY  OF  AMPLITUDE  DISTRIBUTIONS 

We  shall  now  examine  a  number  of  computed  amplitude  probability- 
density  distributions  in  order  to  understand  the  characteristics  impar¬ 
ted  to  such  distributions  by  various  signal  combinations.  We  have  already 
seen  some  distributions  in  the  discussion  of  fading  mechanisms. 

A  single  ray  of  constant  amplitude  will  naturally  produce  an  impulse 
in  amplitude  probability  density  at  the  mean  value  (Fig.  27a).  Combinations 
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Figure  26.  Computed  amplitude  distribution  for  the  random  vector  sum 
of  two  rays. 


(a) 

Amplitudes  =  1.0, 
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(b) 
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(c) 

Amplitudes  =  1.0, 
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of  two  or  three  rays  of  equal  amplitude  produce  the  distributions  shown 
in  Figs.  27b  and  27c.  The  case  of  two  unequal  rays  has  already  been  seen 
in  connection  with  the  discussion  of  polarization  effects  (Fig.  26).  The 
two-ray  distributions  exhibit  a  double  peak,  with  the  spacing  between  the 
peaks  inversely  related  to  the  difference  in  the  strengths  of  the  two 
rays . 

In  Fig.  28,  the  amplitude  distributions  produced  by  three  rays  of 
various  relative  magnitudes  are  plotted.  They  tend  to  exhibit  the  double 
peaked  nature  of  the  two-ray  distributions,  but  are  considerably  smoother, 
and  do  not  have  the  sharp  upper  and  lower  amplitude  limits  of  the  two  ray 
distributions.  The  spacing  of  the  two  peaks  is  seen  to  increase  as  the 
relative  strength  of  the  middle  amplitude  component  increases,  or  as  the 
weakest  component  decreases.  The  clarity  of  the  peaks  is  affected  by  the 
strength  of  the  weakest  component  and  by  the  separation  of  the  peaks. 

As  more  components  of  comparable  magnitude  are  added,  the  character¬ 
istic  features  of  the  distributions  decrease,  approaching  the  Rayleigh 
distribution.  Figure  29a  shows  a  distribution  produced  by  six  rays.  When 
the  number  of  rays  becomes  large,  it  is  not  possible  to  interpret  a  mea¬ 
sured  amplitude  distribution  in  terms  of  known  distributions,  other  than 
to  say  that  a  large  number  of  rays  appears  to  be  present.  An  exception 
occurs  if  the  distribution  is  dominated  by  a  few  strong  rays,  but  other 
rays  of  considerably  lower  amplitude  are  present.  The  weaker  components 
then  may  not  mask  the  features  of  the  distribution,  but  only  smooth  out 
the  features.  Several  distributions  of  this  type  are  shown  in  Fig.  29. 

In  the  discussion  of  polarization  mismatch,  amplitude  distributions 
which  would  be  expected  if  only  one  ordinary  and  one  extraordinary  ray 
should  arrive  at  the  receiver  were  discussed.  If  multiple  ray  paths  ex¬ 
ist,  each  composed  of  an  ordinary  and  an  extraordinary  component,  then 
at  least  four  rays  will  arrive  at  the  receiver.  This  circumstance  will 
increase  the  likelihood  of  having  distributions  which  resemble  the  Ray¬ 
leigh  distribution.  Figures  30a  and  30b  show  examples  of  this  type  of 
distribution  for  equal  ordinary  and  ext rao rd in ary  component  strengths. 

The  representation  used  in  calculating  these  distributions  for  multi-ray 
propagation  with  both  ordinary  and  extraordinary  components  assumes  that 
the  polarization  of  each  of  the  paths  is  independent.  Because  of  the 
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Figure  30.  Computed  amplitude  distributions  for  the  random  vector 
sum  of  rays,  including  the  effect  of  polarization 
mismatch . 

(a)  Ray  polarization-independent,  two  rays,  amplitude  = 
1.0,  0.7,  a2  =  0.232,  cp?  =  o.89. 

(b)  Ray  polarization-independent,  three  rays,  amplitude 
=  1.0,  0.7,  0.4,  a2  =  0.239,  cp?  =  0.97. 

(c)  Ray  polarization-equal,  two  rays,  amplitude  = 

1.0,  0.7,  a2  =  0.453,  y?  =  0.53. 

(d)  Ray  polarization-equal,  three  rays,  amplitude  = 

1.0,  0.7,  0.4,  a 2  =  0.470,  cp?  =  0.71. 
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differences  in  path  lengths  suggested  by  our  model  for  the  rays,  it  is 
anticipated  that  the  resulting  polarizations  would  be  independent.  If 
this  is  not  the  case,  the  polarization  mismatch  for  the  various  paths 
will  be  related.  If  this  assumption  is  used  in  computing  the  expected 
amplitude  distribution,  the  curves  of  Fig.  30c  and  30d  result.  These 
distributions  have  a  very  distinctive  appearance  when  compared  with  the 
other  distributions  we  have  examined,  in  that  they  are  almost  always 
monotonically  decreasing.  Another  characteristic  of  these  distributions 
is  that  the  normalized  variance  usually  exceeds  (4/n)  -1. 

The  amplitude  probability-density  distributions  computed  above  give 
an  idea  of  what  can  be  expected  of  measured  amplitude  distributions  under 
various  conditions.  These  computed  distributions  will  be  compared  with 
experimentally  observed  amplitude  variations  in  the  following  sections. 


D.  RELATION  OF  SIGNAL-STRENGTH  FLUCTUATIONS  TO  PHASE-PATH  FLUCTUATIONS 

In  this  section  a  relation  between  the  occurrence  of  phase-path 
variations  and  fluctuations  in  the  strength  of  the  received  signal  is 
established.  This  relation  is  to  be  expected  on  the  basis  of  the  irreg¬ 
ularity  model  and  its  prediction  of  multipath  interference  during  periods 
of  high  irregularity  activity.  Using  measured  signal-strength  and  phase- 
path  data  acquired  during  the  winter  and  spring  of  1969-1970,  about  ninety 
hours  of  data  were  examined  for  this  purpose.  For  each  1.2  hr  segment  of 
data,  the  rms  phase  fluctuation  s/cp2 

(with  mean  and  linear  trend  removed) 

was  computed.  These  segments  were  grouped  into  ten  intervals  according 
to  the  computed  value  of  For  all  of  the  data  segments  in  each  in¬ 
terval,  the  measured  normalized  variance  a ^  =  [e^/(E  >2]  -  1  of  the 

observed  signal-strength  fluctuation  was  calculated  for  2G0  sec  periods 
(this  period  includes  a  number  of  signal-strength  minima  under  most  con¬ 
ditions).  The  mean  and  the  standard  deviation  of  the  amplitude  variance 
were  computed  for  each  interval  of  *>/cp2.  These  values,  and  the  number  of 
data  segments  in  each  interval,  are  shown  in  Fig.  31.  Although  a  wide 
range  of  signal-strength  fluctuation  is  observed  at  any  measured  value 
of  phase  fluctuation,  a  strong  dependence  between  phase  fluctuations  and 
amplitude  fluctuations  is  seen.  This  indicates  that  the  presence  of  the 
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irregularities  which  produce  the  observed  phase  fluctuations  strong?,y 
influences  the  fading  observed. 

In  order  further  to  verify  the  predicted  effect  of  phase  fluctua¬ 
tions  on  observed  signal  strengths,  a  few  measurements  were  made  over  a 
shorter  path.  It  is  expected  that  for  shorter  paths,  the  effect  of  ir¬ 
regularities  of  a  given  size  and  strength  on  signal-strength  fluctua¬ 
tions  should  not  be  as  great.  This  can  be  seen  from  Eq .  (13),  which 
shows  that  the  expected  number  of  rays  reaching  a  point  on  the  ground 
is  a  function  of  F/A_.  By  decreasing  the  range,  we  reduce  F.  In  ad¬ 
dition,  the  focusing  arising  from  smooth  ionospheric  variations  (Fig.  15) 
will  be  less,  further  reducing  the  effective  value  of  F. 

Simultaneous  observations  were  made  over  the  path  used  for  the  other 
measurements  of  this  work,  and  over  the  path  from  National  Bureau  of  Stan¬ 
dards  Station  WWV  at  Fort  Collins,  Colorado  to  Los  Banos,  California,  in 
order  to  compare  the  influence  of  path  length  on  signal-strength  fluctua¬ 
tions.  Changes  in  phase-path  length  were  also  recorded  simultaneously. 

The  distance  from  Fort  Collins  to  Los  Banos  is  about  1430  km.  The  differ¬ 
ence  in  distance,  plus  ionospheric  focusing,  should  produce  an  effective 
Fresnel-zone  size  for  the  WWV  path  about  half  that  for  the  Arkansas  path. 
This  would  be  expected  to  produce  a  significant  decrease  in  amplitude 
fluctuation  effects. 

Because  no  sounder  was  available  over  the  WWV-Los  Banos  path,  it  was 
not  possible  to  determine  positively  that  only  one  mode  was  propagating 
over  this  path.  However,  ionograms  over  the  Arkansas  path  were  examined 
to  determine  the  most  likely  times  and  frequencies  for  single-mode  propa¬ 
gation  from  Fort  Collins.  This  type  of  propagation  occurred  most  fre¬ 
quently  for  WWV  at  15  MHz.  A  single-mode  frequency  for  the  Arkansas  path 
was  chosen  to  provide  a  reflection  height  which  was  comparable  to  that  on 
the  WWV  path  based  on  nominal  ionospheric  conditions.  This  was  usually 
about  24  MHz.  (The  difference  in  frequencies  will  cause  the  Fresnel-zone 
length  for  WWV  to  be  slightly  larger  than  halx  the  Fresnel-zone  length 
f  or  Bearden . ) 

In  all  observations,  fading  on  the  WWV  path  was  significantly  less 
than  on  the  Arkansas  path.  At  times  the  WWV  path  exhibited  almost  no 
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noticeable  fading.  Figure  32  is  typical  of  the  data  obtained  under  such 
conditions.  The  phase-path  variations  are  similar  in  magnitude,  although 
uncorrelated.  The  signal  fluctuations  on  the  W\W  path  are  considerably 
weaker.  In  addition,  the  measured  amplitude  density  distributions  for 
the  period  of  the  records  (Fig.  33)  shows  that  the  WWV  signal  is  much 
more  constant  in  amplitude.  Again,  this  verifies  our  interpretation  of 
the  phase  fluctuations  observed. 


E.  MEASURED  AMPLITUDE  HISTOGRAMS 

We  shall  now  examine  a  number  of  amplitude  histograms  calculated 
from  the  measured  signal-strength  records,  and  attempt  to  infer  some 
characteristics  of  the  rays  producing  the  observed  distributions.  A 
major  difficulty  in  this  task  is  presented  by  the  discreteness  of  the 
irregularities  and  the  small  number  of  rays.  Under  these  conditions, 
the  number  and  amplitudes  of  the  rays  are  expected  to  vary  significantly 
with  time.  For  irregularities  having  a  20  min  period,  it  is  net  likely 
that  the  situation  will  remain  stationary  for  more  than  a  few  minutes. 

For  observations  to  represent  conditions  at  any  one  particular  time,  the 
conditions  should  remain  relatively  constant  for  the  period  of  the  mea¬ 
surement.  The  measuring  period  will  therefore  be  a  compromise  between  a 
period  short  enough  to  provide  stationarity  and  a  period  long  enough  to 
permit  obtaining  the  reasonably  large  number  of  samples  required  to  pro¬ 
duce  valid  distributions. 

Measurements  of  amplitude  were  made  at  intervals  of  0.1  to  1.0  sec. 

There  is  a  limit  to  the  amount  of  information  to  be  gained  by  sampling 

more  rapidly.  If  the  observed  distribution  is  to  be  valid,  the  rays  must 

have  doppler  shifts  sufficiently  large  and  sufficiently  different  from 

each  other  that  a  random  phasor  sum  of  the  rays  will  be  observed.  The 

effects  of  different  sampling  intervals  will  be  discussed  together  with 

2 

the  observed  amplitude  distributions.  Also,  the  variance  a  of  the 

observed  amplitude  fluctuations  will  be  noted  for  comparison  with  those 

2 

predicted  in  the  computed  distributions.  The  values  of  a  are  indica¬ 
ted  in  the  captions  for  each  of  the  computed  and  measured  distributions. 
If  the  distributions  result  from  a  random  phasor  sum  of  sine  waves  whose 
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Figure  32.  Observed  s ignal -s trength  and  phase-path-length 
fluctuations  over  two  paths  at  the  same  time. 

(u)  Signal  strength,  WWV-15,  Fort  Collins, 

Colorado  to  Los  Banos,  Calif.  (1430  km). 

(c)  Phase-path-length  f luctuations--For t  Collins. 

Continued  on  next  page. 
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Figure  32.  Continued. 


(b)  Signal  strength--22 . 81  MHz,  Bearden, 

Arkansas  to  Los  Banos,  Calif.  (2600  km) 
(d)  Phase-path-length  f luctuations--Bearden . 
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Figure  33.  Observed  amplitude  density  distributions  W(E/E) 
for  the  period  shown  in  the  records  of  Fig.  33. 

o 

(a)  Fort  Collins  path.  -  a  =  0.081  . 

2 

(b)  Bearden,  Arkansas  path.  -  a  =  0.204  . 
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amplitudes  do  not  vary,  the  variance  a  is  restricted  to  lie  in  the 
range  from  0  to  0.273  [(4/n)  -1],  Greater  values  may  occur  if  the  sig¬ 

nal  energy  changes  during  the  analysis  interval.  Such  variations  in  en¬ 
ergy  would  be  expected,  based  on  our  irregularity  model. 

The  examination  of  a  large  number  of  measured  amplitude  distributions 
reveals  some  general  characteristics.  First,  no  distribution  curve  shape 
may  be  said  to  be  typical.  The  shape  of  the  distributions  is  a  strong 
function  of  the  analysis  interval  and  the  nature  of  the  signal  fluctua¬ 
tions  during  the  analysis  interval.  The  second  general  characteristic 

2 

of  the  data  is  that  a  frequently  exceeds  0.273,  indicating  that  for 
analysis  times  long  enough  to  produce  valid  random  phasor  sums,  signifi¬ 
cant  energy  changes  occur.  Finally,  only  when  fading  is  quite  rapid  and 
severe  do  the  measured  distributions  resemble  any  of  the  Rice  distribu¬ 
tions  shown  in  Figs.  21  and  23,  and  then  only  the  p  =  0  distribution, 
which  is  equivalent  to  the  computed  distributions  for  more  than  a  few 
rays  of  comparable  magnitude. 

A  frequently  noted  characteristic  of  the  observed  amplitude  distri¬ 
butions  is  a  double-peaked  form.  Three  such  distributions  are  shown  in 
Fig.  34,  which  displays  the  observed  amplitude  density,  W(E/E),  versus 
E/E.  These  distributions  are  measured  from  data  taken  during  10  to  15 
min  intervals .  These  double-peaked  curves  resemble  the  2-ray  and  3-ray 
distributions  of  Figs.  26  and  28.  They  differ  from  the  latter  distribu¬ 
tions  in  having  more  elongated  high-amplitude  "tails".  Such  long  tails 
are  frequently  observed,  and  are  assumed  to  be  caused  by  changes  in  the 
amplitudes  of  the  rays  being  received. 

In  an  effort  to  verify  this  assumption,  some  data  were  analyzed  for 
a  2.5  min  period.  These  data  were  sampled  at  10  samples  per  sec,  so  that 
enough  points  to  produce  a  smooth  distribution  could  be  obtained.  Records 
such  as  shown  in  Fig.  35a  were  produced  under  these  conditions.  This 
distribution  has  a  sharp  upper  limit.  It  does  not,  however,  appear  to 
be  a  satisfactory  random  vector  sum,  and  while  ray  amplitudes  are  prob¬ 
ably  relatively  constant  for  this  interval  of  2.5  min,  the  frequency 
shifts  of  the  various  rays  are  apparently  not  sufficient  to  produce  the 
desired  random  vector  sum.  Figure  35b  shows  the  distribution  produced 
by  analyzing  a  5  min  period  which  included  the  period  used  for  Fig.  35a. 
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Figure  35.  Observed  amplitude  density  distributions 
for  short  analysis  times. 

2 

(a)  Analysis  time  =  2,5  min.  -  a  =  0.194 

2 

(b)  Analysis  time  =  5-.0  min.  -  a  =  0.238 


Note  that  the  tail  of  the  distribution  is  growing,  indicating  a  change 
in  the  amplitudes  of  some  of  the  component  rays. 

To  further  investigate  the  effect  of  changing  ray  amplitudes,  a 
distribution  was  computed  for  which  one  ray  was  present  for  half  of  the 
samples  and  three  rays  were  present  for  the  other  half.  Such  a  stepwise 
change  is  not  realistic,  but  the  resulting  computed  distribution,  shown 
in  Fig.  36a,  gives  an  insight  into  the  effect  of  changing  energy.  A 
measured  distribution  exhibiting  a  similar  effect  is  shown  in  Fig.  36b. 

All  of  the  above  measured  distributions  were  taken  during  periods 
of  relatively  rapid  fading,  as  observed  on  the  signal-strength  records. 

If  the  observed  fading  were  the  result  of  interference  between  two  rays 
of  about  equal  amplitude,  such  as  would  be  the  expected  result  of  magne¬ 
toionic  splitting,  distributions  similar  to  those  of  Figs.  26,  30a  and 
30b  should  be  observed.  Records  such  as  these,  however,  are  not  observed 
during  periods  of  rapid  signal-strength  fluctuation,  indicating  that 
Faraday  rotation  is  not  the  dominant  mechanism  in  producing  these  fades. 
Observations  made  over  long  periods  of  time  when  little  fading  is  present, 
do  however  appear  quite  similar  to  the  computed  distributions  for  two 
rays  of  similar  strength.  Figures  37a  and  37b  show  two  such  measured 
distributions,  where  the  analysis  time  is  20  min.  This  indicates  that 
the  effects  of  Faraday  rotation  are  seen  over  longer  time  intervals,  but 
not  for  short  intervals.  This  is  consistent  with  predicted  behavior  for 
the  ordinary  and  extraordinary  rays. 

Additional  examples  of  some  observed  distributions  are  shown  in  Fig. 
38.  These  distributions  were  measured  during  three  consecutive  10  min 

periods  under  conditions  of  moderate  fading.  They  show  the  change  in 

2 

the  dominant  rays  during  the  interval.  In  all  cases  a  is  less  than 
0.273,  indicating  that  energy  was  relatively  constant  during  each  10  min 
period.  Figure  38d  shows  the  amplitude  distribution  for  the  entire  pe¬ 
riod,  indicating  the  loss  of  information  due  to  changes  in  the  ray  struc- 

2 

ture.  Also,  for  this  record,  a  is  greater  than  0.273, 

Finally,  we  show  three  more  measured  distributions  in  Fig.  39.  These 
are  the  result  of  analyzing  20  min  segments  of  data  during  a  one  and  a 
half  hour  period.  The  signal-strength  record  for  this  period  is  shown 
in  Fig.  39d .  The  three  distributions  correspond  to  periods  of  time  when 


COMPUTED  DENSITY  4>( 


Figure  36.  Amplitude  distributions  exhibiting  the  effect  of  changing 
signal  energy. 

(a)  Computed  distribution  with  one  ray  (amplitude  =  1.0) 

present  for  half  of  the  analysis  interval  and  three 
rays  (amplitude  -  1.0,  0.85,  0.7)  present  for  the 
other  half  of  the  analysis  interval.  -  -  0.288, 

9  =  0.79  . 

2 

(b)  Measured  distribution.  -  a  =  0.26  . 
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Figure  37.  Measured  amplitude  distributions  showing 
the  effects  of  Faraday  rotation. 

2 

(a)  Analysis  time  =  20  min.  -  a  =0.14 

2 

(b)  Analysis  time  =  20  min.  -  a  =  0.20 


Figure  38.  Measured  amplitude  distribution  for  three  consecutive 
10-min  intervals,  and  for  the  whole  30-min  interval, 
showing  the  effect  of  analysis  time  on  observed 


>.  -S'  w/f . 


the  signal-strength  record  had  different  appearances.  Figure  39a  appears 
to  resemble  the  amplitude  distribution  produced  by  two  or  three  rays  (Fig 
26).  During  the  corresponding  interval  the  signal-strength  record  shows 
a  moderate  amount  of  fading.  Figure  39b,  the  distribution  for  a  more  dis 
turbed  period  on  the  signal  strength  record,  more  closely  resembles  compu 
ted  distributions  for  a  larger  number  of  rays.  Figure  39c  shows  a  very 
strong  component  indicating  that  a  single  ray  was  prominent  during  this 
period.  The  signal  strength  during  this  latter  time  was  very  stable. 


F.  CONCLUSIONS  FROM  SIGNAL-STRENGTH  ANALYSIS 

The  observed  amplitude  distribution  data  presented  in  the  above  ex¬ 
amples  are  similar  to  a  majority  of  the  observed  amplitude  distributions 
examined.  Analysis  times  were  chosen  to  emphasize  certain  features  of 
the  data.  For  all  data  taken  during  the  course  of  the  experiments,  sim¬ 
ilar  distributions  were  produced  when  analysis  intervals  were  chosen  long 
enough  to  give  reasonably  smooth  distributions,  and  short  enough  to  avoid 
encompassing  widely  varying  situations  which  produce  distributions  ap¬ 
proaching  the  Rayleigh.  From  these  analyses,  the  following  conclusions 
are  drawn : 


1.  Signal-strength  fluctuations  observed  for  single-mode  iono¬ 
spheric  propagation  under  quiet  ionospheric  conditions  are 
adequately  described  in  terms  of  the  multipath  interference 
resulting  from  vertical  ray  bending  caused  by  large,  weak 
ionospheric  irregularities. 

2.  Under  most  conditions,  fading  due  to  0-X  interference  is 
not  significant  compared  with  that  arising  from  the  vertical 
ray  bending  caused  by  the  irregularities  being  investigated. 

3.  Measured  amplitude  histograms  for  ionospherically  propagated 
rays  are  highly  time-dependent  because  of  the  variability  of 
the  causative  mechanism. 

4.  Weak  scattering  caused  by  a  large  number  of  small  scatterers 
does  not  produce  the  major  portion  of  the  fading  observed  on 
single-mode  one-hop  ionospheric  propagation.  We  do  not  know 
whether  such  weak  scattering  occurs  in  some  smaller  degree 
which  is  masked  by  the  multipath-interference  effects.  Ex¬ 
periments  to  answer  this  question  must  consider  the  effects 


of  multipath  such  as  we  observe.  Such  experiments  would  promise 
most  success  if  performed  for  shorter  paths,  where  multi-ray  in¬ 
terference  is  less,  and  during  periods  of  low  fluctuation  in 
phase-path  length. 
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Figure  40.  Doppler  correlation  versus  time  lag 

for  wide-spaced  phase-path  measurements 
made  with  receiving  antennas  at  Los 
Banos,  Calif,  and  on  an  airplane. 
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V .  OTHER  CHARACTERISTICS  OF  THE  IRREGULARITIES 

A.  HORIZONTAL  MOTION  OF  IRREGULARITIES 

Some  of  the  phase-path  data  acquired  in  this  work  have  been  examined 
with  the  object  of  determining  whether  the  observed  irregularities  exhibit 
any  preferred  direction  and  velocity  of  motion.  The  observational  geometry 
used  in  this  work  was  sensitive  to  horizontal  motions  in  the  north-south 
direction,  and  also  to  vertical  motions  (when  more  than  one  radio  frequency 
was  used  simultaneously). 

The  records  obtained  from  the  wide-spaced  CW  measurements,  using  re¬ 
ceivers  at  Los  Banos  and  on  an  airplane,  were  examined  in  an  attempt  to 
identify  a  consistent  horizontal  velocity  component  of  the  irregularities. 
The  aircraft  doppler  measurements  taken  at  spacings  of  60,  65  and  70  km 
were  correlated  [as  in  Eqs .  (7)  and  (8)]  with  the  fixed  receiver-site 
measurements  using  time  lags  of  ±50  sec  to  ±200  sec  for  the  fixed  site, 
to  determine  whether  the  highest  correlation  occurred  at  a  time  lag  other 
than  zero.  The  results  of  these  correlations  are  shown  in  Fig.  40.  In 
all  cases,  the  correlation  decreased  for  time  lags  other  than  zero,  indi¬ 
cating  no  preferred  horizontal  motion. 

The  above  test  involved  averaging  data  over  eight  observations,  and 
therefore  would  sense  only  a  preferred  direction  of  motion.  It  is  possi¬ 
ble  that  each  of  the  irregularities  is  moving,  but  that  the  directions 
of  motion  vary  widely.  In  order  to  observe  the  motion  of  individual  ir¬ 
regularities,  use  was  made  of  the  data  gathered  during  the  winter  and 
spring  of  1969-1970.  Charges  in  phase-path  length  measured  on  two  ver¬ 
tical  monopoles  spaced  2.5  km  apart  were  cross-correlated  to  determine 
the  time  shift  of  maximum  correlation.  Figure  41a  shows  a  typical  result 
of  this  correlation.  The  phase  of  the  signals  at  this  spacing  was  highly 
correlated  (greater  than  0.97  in  all  cases).  The  maximum  phase-path  cor¬ 
relation  pp(At)  occurs  for  a  time  shift  of  0  ±  5  sec  for  all  records 
examined.  Because  of  the  slow  decorrelation  of  phase  with  time,  it  is 
difficult  to  assign  a  value  of  north-south  velocity  component  from  these 
records,  and  no  attempt  was  made  to  do  so. 

Similar  correlations  were  performed  on  the  amplitude  data  available 
from  the  experiments.  As  shown  in  Fig.  41b,  the  amplitude  correlation 
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Figure  41.  Cross  correlation  of  signals  on  antennas  spaced  2.5  km  apart, 
(a)  Phase-path  correlation,  Pp(At)  . 

Continued  on  next  page. 


(b)  Envelope  amplitude  corre 
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Ef  PE<.J\t)  decreases  much  more  rapidly  with  time  than  does  the  phase-path 

*3 

1 

( 

correlation.  Offsets  of  the  peak  amplitude  correlation  (from  zero  lag) 

i  ranged  from  0  sec  to  z20  sec.  The  interpretation  of  these  time  lags  is 

■ 

1 

1 

\  difficult  for  two  reasons.  First,  only  one  component  of  the  velocity  is 

1 

observed.  This  allows  us  only  to  place  a  lower  limit  on  the  velocity  of 

■ 

the  fading  pattern.  Second,  horizontal  motion  of  the  fading  pattern  can 

3 

i; 

result  from  vertical  motion  and  deformation  of  the  irregularities  as  well 

fj 

as  from  their  horizontal  motion. 

, 

If  we  assume  that  horizontal  irregularity  motion,  rather  than  verti- 

cal  motion  or  deformation,  is  responsible  for  the  motion  of  the  fading 

pattern,  we  find  (based  on  the  antenna  spacing  and  time  lag  given  above) 

a  lower  limit  for  the  horizontal  velocity  of  the  irregularity  of  about 

65  m  per  sec . 

\ 

:  B.  VERTICAL  MOTION  OF  IRREGULARITIES 

! 

p 

By  making  use  of  simultaneous  transmissions  on  two  frequencies  over 

’i 

tj 

j  the  Arkansas-Los  Banos  path,  changes  in  phase-path  length  could  be  mea- 

;  surod  on  two  vertically  separated  paths  simultaneously.  Approximately 

thirty  hours  of  such  two-frequency  data  taken  during  the  winter-spring 

1 

1969-1970  period  were  analyzed  to  determine  whether  the  observed  irregu- 

['  larities  moved  vertically.  The  analysis  was  performed  by  cross-correla- 

E 

ting  the  phase-path  variation  records  for  the  two  paths.  Data  segments 

■ 

of  1.2  hours  were  used  for  each  correlation.  Frequency  separations  be- 

i 

*  tween  100  kHz  and  5.1  MHz  were  used.  A  typical  record  for  a  2  MHz  fre- 

1 

1  quency  spacing  is  shown  in  Fig.  42.  The  peak  correlation  occurs  at  a 

time  lag  of  240  sec  for  the  lower-frequency  path.  Assuming  the  vei’tical 

i 

ray  spacings  given  in  Fig.  8,  this  corresponds  to  a  downward  velocity  of 

1 

about  33  m/sec.  All  of  the  data  analyzed  in  this  manner  (shown  in  Fig. 

B 

43)  indicated  a  downward  velocity  of  between  30  and  40  m/sec .  We  there- 

B 

fore  have  considerable  confidence  in  this  measurement. 
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DIURNAL  AND  SEASONAL  VARIATIONS 


The  measuring  technique  used  in  the  course  of  this  work  could  be 
regularly  employed  only  during  winter  daytime  periods  for  the  path  in¬ 
volved.  To  determine  any  diurnal  variation  which  might  be  observable 
during  these  limited  periods  (about  10:00  am  to  4:00  pm  local  time  at 
midpath)  the  rms  phase  fluctuations  observed  for  each  1.2  hr  segment 
of  data  analyzed  in  Chapter  III,  Section  D  (Fig.  31)  are  grouped  into 
1  hour  intervals  (according  to  the  time  of  day  at  the  middle  of  the 
segment)  and  an  average  is  calculated.  These  averages  and  their  stan¬ 
dard  deviations  are  shown  in  Fig.  44.  They  indicate  that  large  phase 
fluctuations  are  most  common  in  the  morning,  and  least  common  near  noon. 

To  examine  the  seasonal  variations  of  \|cj?  this  quantity  is  aver¬ 
aged  for  each  day  of  available  data.  These  averages  are  shown  in  Fig. 
45,  for  the  period  from  November  1969  through  April  1970.  Single-mode 
phase-path  measurements  are  not  possible  on  a  regular  basis  over  the 
Bearden-Los  Banos  path  at  other  times  of  the  year.  For  the  period  of 
the  data,  a  small  decrease  in  the  magnitude  of  the  phase-path-length 
fluctuations  is  observed  during  the  latter  portion  of  the  period. 


D .  ANTENNA  PERFORMANCE 


If  it  is  desired  to  make  use  of  very  large  apertures  (greater  than 
one  km)  for  HF  communication  systems,  in  order  to  increase  antenna  gain 
and  reduce  interference,  the  irregularities  described  in  this  work  will 
strongly  influence  the  design  and  performance  of  such  antennas.  The 
strength  and  horizontal  scale  size  of  the  irregularities  will  determine 
the  wave-front  curvature  for  rays  arriving  at  the  receiving  antenna. 

This  will  limit  the  length  of  antennas  which  may  be  employed  without 
compensation  for  the  varying  curvature  of  the  wave  fronts.  It  is  found, 
however,  that  antennas  having  lengths  up  to  about  20  km  will  be  useful 
under  many  circumstances  without  compensation  for  wave-front  curvature. 
Such  antennas  must,  however,  be  able  to  track  the  changing  apparent 
azimuthal  angle  of  arrival  of  incoming  rays. 
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Figure  44.  Variation  of  phase  fluctuations  with  time 
of  day. 
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Figure  45.  Variation  of  phase  fluctuations  with  time 
of  year. 
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S Lgua 1 -s t rength  fluctuations  will  produce  fading  on  any  antenna 
that  is  short  compared  with  the  dimensions  of  the  ray  interference 
pattern  on  the  ground .  Measurements  made  during  the  course  of  this 
work  show  that  the  correlation  coefficient  for  amplitude  fluctuations 
on  two  spaced  antennas  perpendicular  to  the  ray  path  is  0.8  to  0.95  at 
spactngs  of  360  m  and  0.4  to  0.8  at  spacings  of  2.5  km  (Fig.  41b). 
larger  horizontal  spacings,  vertical  spacings,  or  spacings  which  are 
not  normal  to  the  ray  path  may  provide  more  rapid  decreases  in  the 
amplitude  correlation  with  distance.  If  the  antenna  can  be  made  large 
enough  so  that  it  is  composed  of  a  number  of  cells  (subdivisions)  each 
of  which  observes  amplitudes  independent  of  those  observed  by  the  other 
cells,  the  antenna  performance  may  be  predicted  using  the  methods  of 
Dracewell  (’60).  With  this  method,  the  field  at  the  receiving  antenna, 
E,  is  represented  as  the  sum'  of  the  strongest  ray  arriving  at  the  an¬ 


tenna,  E^,  and  a  superimposed  random  component,  where  (•  is 


the  complex  deviation  factor  of  the  field  from  its  average  value.  The 
efficiency  of  the  antenna  is  then  given  by 


S  = 


1  +  var  £ 


(20) 


In  Chapter  III  the  variance  of  the  amplitude  fluctuations  (a  ) 
for  each  measured  and  computed  distribution,  and  the  variance  of  the 
phase  fluctuations  (cp2)  for  the  computed  distributions  were  calcula¬ 
ted.  These  variances  may  be  used  to  predict  the  efficiency  of  large 
antenna  systems  if  the  antenna  is  large  enough  that  the  required  inde¬ 
pendent  cells  are  present.  Fading  would  not  be  observed  with  such  an 
antenna,  although  its  efficiency  would  be  reduced  by  the  multiple  in¬ 
cident  rays. 


E . 


BANDWIDTH  LIMITATIONS 


The  multiray  representation  of  the  field  resulting  from  propaga¬ 
tion  through  ionospheric  irregularities  may  also  be  used  in  describing 
the  bandwidth  characteristics  of  the  channel.  At  any  given  frequency, 
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the  channel  transfer  function  is  described  by  the  relative  amplitudes  of 
the  rays  and  by  their  differential  time  delay,  AP/c .  The  manner  in  which 
these  quantities  vary  with  frequency  has  not  been  investigated  in  this 
work . 

If  the  directivity  of  an  antenna  can  be  made  high  enough  so  that 
only  one  ray  is  accepted,  the  interfering  rays  will  not  be  admitted  to  > 

\ 

the  receiver.  This  would  eliminate  the  fading  and  the  bandwidth  limita-  \ 
tions  described  above.  In  the  vertical  (elevation)  plane,  rays  will 
originate  at  heights  which  differ  by  about  15  km  at  the  reflection  height. 
For  a  path  of  the  length  used  in  our  measurements,  these  rays  will  arrive 
at  the  receiving  antenna  with  a  separation  in  elevation  angle  of  about  a 
half  degree.  A  vertical  aperture  of  over  1  km  would  be  required  to  sepa¬ 
rate  these  rays.  Therefore,  it  does  not  appear  practical  to  reject  inter¬ 
fering  rays  by  the  use  of  large  vertical  apertures.  However,  large  hori¬ 
zontal  antennas  might  be  able  to  reject  multipath  rays,  if  the  rays  have 
sufficiently  different  azimuthal  angles  of  arrival. 
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VI.  SUMMARY  AND  RECOMMENDATIONS 


A .  SUMMARY 

High-frequency  ionospheric  radio  propagation  over  a  one-hop  lower- 
ray  F-layer  path  has  been  studied  under  "quiet"  ionospheric  conditions, 

i.e.,  during  periods  when  high-resolution  oblique-incidence  soundings 
over  the  path  showed  no  evidence  of  ionospheric  irregularities.  Changes 
in  the  phase-path  length  of  rays  propagating  over  the  path  were  measured 
using  a  technique  which  is  sensitive  to  changes  occurring  as  slowly  as 
ten  meters  per  hour.  Under  quiet  conditions,  the  technique  allows  mea¬ 
suring  ionospheric  refractive-index  irregularities  with  a  sensitivity 
about  an  order  of  magnitude  better  than  that  of  other  methods,  e.g., 
pulse  sounders,  satellite  Faraday  rotation  measurements  and  doppler 
spectrum  analysis. 

Analysis  of  the  data  produced  in  these  measurements  revealed  a 
class  of  weak  s lowly -vary ing  irregularities  which  were  present  during 
all  of  the  observations.  The  data  were  taken  during  daylight  hours  in 
the  spring  of  1969  and  the  winter  and  spring  of  1969-1970,  over  a  2600 
km  mid-latitude  east-west  path.  A  variety  of  experimental  configura¬ 
tions  were  employed.  The  major  physical  characteristics  of  the  irre¬ 
gularities  are: 

1.  Horizontal  scale  size  of  about  35  km. 

2.  Vertical  scole  size  of  about  15  km. 

3.  Electron-density  changes  of  0.6  to  6.0  per  cent. 

4.  Horizontal  velocity  of  65  m/sec  (lower  limit) — no  preferred 
direction . 

5.  Vertical  phase  velocity  of  35  m/ sec --downward . 

Using  these  parameters,  we  may  classify  the  new  irregularities  as 
smaller,  weaker,  and  somewhat  slower-moving  than  those  studied  by  Munro, 
Titheridge  and  others.  Their  continuous  presence  and  comparatively  large 
effects  make  them  quite  significant  in  ionospheric  radio  propagation  over 
long  one-hop  F-layer  paths. 


A  model  was  developed,  and  was  shown  to  predict  correctly  the 
amplitude  fluctuations  most  commonly  observed  for  long-distance  sin¬ 
gle-mode  HF  propagation.  These  fluctuations  had  not  been  satisfacto¬ 
rily  explained  previously.  The  irregularities  also  impose  fundamental 
limits  on  the  spatial  resolution  and  bandwidth  which  may  be  achievable 
in  HF  ionospheric  propagation. 

No  attempt  has  been  made  to  determine  the  origin  of  these  irregu¬ 
larities,  but  their  characteristics  are  consistent  with  effects  which 
might  be  produced  by  internal  gravity  waves.  Georges  (’67)  in  studying 
large-  and  medium-scale  traveling  ionospheric  disturbances  points  out 
that  preferred  periods  for  gravity  waves  in  the  F-region  are  close  to 
20  minutes.  He  also  shows  that  for  internal  gravity-wave  energy  propa¬ 
gating  vertically  upward,  the  phase  velocity  of  the  resulting  distur¬ 
bances  is  vertically  downward.  If  the  irregularities  observed  in  this 
investigation  result  from  internal  gravity  waves,  there  is  reason  to 
believe  that  the  source  of  energy  producing  them  lies  below  the  F  region, 
and  possibly  within  the  upper  atmosphere.  Georges  notes  that  thunder¬ 
cloud  tops  are  observed  to  oscillate  with  periods  of  about  20  min,  and 
that  attempts  have  been  made  by  others  to  detect  ionospheric  effects 
which  might  be  produced  by  these  oscillations.  The  phase-path  measure¬ 
ment  method  used  in  this  work  offers  considerable  potential  for  observ¬ 
ing  such  a  relation  if  it  exists. 

B.  RECOMMENDATIONS 

The  bandwidth  limitations  imposed  by  weak  ionospheric  irregularities 
are  not  yet  fully  understood.  They  are,  however,  an  important  result  of 
propagation  through  such  irregularities.  Additional  information  concern¬ 
ing  the  characteristics  and  occurrence  of  these  irregularities  is  needed 
in  order  to  form  a  complete  description  of  their  effect  on  HF  propagation 
under  varying  circumstances.  The  stable-frequency  phase-path  measuring 
system  employed  in  this  work  would  be  particularly  suitable  for  additional 
investigations  of  the  motion,  the  solar-cycle  dependence,  and  the  seasonal 
diurnal  and  global  frequency  of  occurrence  of  these  irregularities.  For 
use  in  potential  development  of  large  antenna  systems,  more  detailed  infor 
raation  on  the  azimuthal-angle  distribution  of  multiple  rays  is  needed. 


by  these  irregularities  is  multipath  interference  within  a  single  mode. 
Since  the  extent  of  this  interference  is  a  function  of  the  distance  from 
the  scattering  region  to  the  receiver,  it  may  be  possible  to  optimize  the 
design  of  a  high-performance  communication  system  by  selecting  paths  suf¬ 
ficiently  short  that  such  interference  is  minimized  to  the  greatest  possi¬ 
ble  extent.  More  data  is  required  on  the  occurrence  and  global  distribu¬ 
tion  of  the  irregularities  if  such  an  optimization  is  to  be  used  in  prac¬ 
tice  . 
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APPENDICES 


APPENDIX  A 


CW  PHASE -PATH  MEASUREMENTS 

A.  THE  BASIC  MEASUREMENT 

The  phase-path  measurements  of  this  work  were  made  using  CW  trans¬ 
missions  from  stable-frequency  transmitters.  The  transmitted  frequencies 
were  synthesized  from  a  cesium-beam  frequency  standard  (or  in  some  cases, 
from  a  highly  stable  crystal  frequency  standard).  The  signals  were  re¬ 
ceived  using  standard  HF  receivers.  All  of  the  conversion  injection  fre¬ 
quencies  of  these  receivers  were  synthesized  from  cesium-beam  frequency 
standards.  The  audio-frequency  output  of  the  receivers  was  sampled  by  a 
computer  with  an  analog-to-digital  converter.  Zero  crossings  of  the  au¬ 
dio  signal  were  counted  continuously  from  the  start  of  each  data-gather- 
ing  run.  At  each  "sampling  interval"  the  computer  subtracted  a  number, 
n,  from  the  running  total  of  zero  crossings.  The  value  of  n  was 

n  =  2f  t  , 
a 

where  f  is  the  nominal  audio  output  frequency  of  the  receiver,  as 

cl 

determined  by  the  receiver  conversion  frequencies  used,  and  t  is  the 
sampling  interval  used.  Sampling  intervals  between  1/10  and  1  sec  were 
used.  The  sampling  rate  and  the  length  of  the  sampling  interval  were 
controlled  by  pulses  derived  from  digital  dividers  clocked  by  the  cesium- 
beam  frequency  standard. 

In  addition,  at  the  end  of  each  sampling  interval  the  instantaneous 
phase  of  the  audio  signal  was  calculated  by  measuring  the  fractional  part 
of  a  half -cycle  as 


f  ract . 


where  t^  is  the  time  between  the  end  of  the  interval  and  the  zero  cross¬ 
ing  immediately  prior  to  the  end  of  the  sampling  interval,  and  t ^  is 
the  time  between  the  zero  crossings  immediately  preceding  and  immediately 
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following  the  end  oi'  the  sampling  interval.  The  sum  of  the  running  total 
of  zero  crossings  and  the  measured  fraction  is  the  desired  measurement  of 
phase-path  length  in  units  of  A/2,  where  A  is  the  wavelength  of  the 
transmitted  frequency.  Using  this  technique,  errors  in  the  measurement 
of  the  fractional  part  of  a  half-cycle  do  not  accumulate,  as  the  output 
at  each  sampling  interval  is  the  running  total  of  zero  crossings  for  the 
entire  measurement  period  (an  integer  number)  plus  the  fraction  measured 
at  the  end  of  the  specific  interval. 

The  envelope  amplitude  of  the  signal  was  determined  by  measuring  the 
slope  of  the  audio  signal  at  the  zero  crossings  immediately  preceding  and 
following  the  end  of  the  sampling  interval.  For  narrow-band  signals,  this 
slope  is  directly  proportional  to  the  envelope  amplitude  of  the  signal. 

Errors  in  these  measurements  are  introduced  by  instabilities  in  the 
frequency  standards  and  in  the  other  equipment  (e.g.,  receiver  phase  drift, 
digital-circuit  jitter).  As  stated  in  the  body  of  the  report,  errors  in 
the  frequency  standards  introduce  an  error  of  about  10  m/hr  in  the  mea¬ 
surement  of  phase-path  length.  Tests  were  made  to  determine  the  total 
system  error  by  introducing  synthesized  standard  frequencies  directly 
into  the  input  of  the  receivers.  The  long-term  phase  drift  observed  was 
less  than  10  m/hr  (equivalent  phase-path  length)  when  the  receiver  input 
was  derived  from  the  same  frequency  standard  used  to  synthesize  the  re¬ 
ceiver  injection  frequencies,  or  from  another  independent  cesium-beam 
frequency  standard.  When  a  stable  crystal  frequency  standard  was  used 
to  supply  the  receiver  input,  a  linear  long-term  drift  was  observed  due 
to  the  frequency  error  of  the  crystal  standard.  Departures  from  a  linear 
drift  were  less  than  20  m/hr. 

Short-term  phase  errors,  arising  primarily  from  noise  and  digital 
circuit  jitter,  were  less  than  2  deg  at  the  maximum  receiver  input  level, 
and  less  than  10  deg  (~  1/3  m)  at  a  receiver  input  40  dB  below  the  maxi¬ 
mum  level.  The  amplitude  measurement  was  accurate  to  within  1  dB  over 
the  40  dB  range  of  input  levels. 

During  the  first  series  of  measurements,  involving  the  airborne  and 
Los  Banos  receiver  sites,  some  additional  short-term  errors  were  intro¬ 
duced  by  the  intermediate  analog  magnetic  tape  recording.  These  errors 
are  discussed  in  the  next  section. 
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B.  WIDE-SPACED  MEASUREMENTS 


Initial  determinations  of  the  strength  and  horizontal  size  of  the 
observed  weak  irregularities  were  made  using  a  transmitter  in  Bearden, 
Arkansas,  a  fixed  receiving  site  at  Los  Banos,  California,  and  a  receiv¬ 
ing  system  in  an  aircraft.  The  path  length  is  about  2600  km.  This  con¬ 
figuration  was  used  to  allow  a  relatively  rapid  determination  of  the  ma¬ 
jor  characteristics  of  the  received  wavefronts.  Flights  were  made  in  a 
generally  north-south  line,  from  30  to  170  km  north  of  the  fixed  receiv¬ 
ing  site  (see  Fig.  46).  Each  flight  took  between  20  and  40  min,  and 
covered  from  100  to  140  km.  Ten  flights  were  made,  of  which  the  eight 
most  successful  were  used  to  determine  the  parameters  of  the  irregular¬ 
ities.  This  section  describes  the  equipment  and  methods  used  in  making 
these  measurements . 

A  DC-3  aircraft  was  used  for  the  experiment.  It  was  equipped  with 
two  monopole  antennas,  approximately  1/4  wavelength  long.  These  were 
tuned  so  that  a  low  VSWR  was  obtained  over  frequencies  from  20  to  26  MHz. 
The  same  signal  was  received  on  both  of  the  antennas  and  separately  re¬ 
corded.  The  x’esults  of  processing  data  from  both  antennas  were  compared 
to  eliminate  any  ambiguities  resulting  from  fades  due  to  ground  reflec¬ 
tions  or  incoming  polarization  changes.  Generally,  the  data  from  the 
vertical  antenna  were  used.  Figure  47  is  a  photograph  of  the  exterior 
of  the  aircraft.  The  antennas  may  be  seen  at  the  rear. 

Figure  4P/  shows  the  interior  of  the  aircraft  with  the  equipment  in¬ 
stalled.  All  receivers  were  modified  so  that  all  conversion  frequencies 
were  synthesized  from  the  cesium  standard.  The  receiver  audio  outputs 
were  recorded  on  analog  magnetic  tape.  In  addition,  a  50  kHz  signal  de¬ 
rived  from  the  cesium  standard  was  recorded  on  the  tape.  This  signal  was 
used  during  data  reduction  to  control  the  sampling  of  the  data  channels, 
thereby  greatly  reducing  the  effects  of  short-term  tape  speed  fluctua¬ 
tions,  and  eliminating  the  effects  of  long-term  fluctuations.  A  one 
pulse  per  second  (pps)  signal  derived  from  the  cesium  standard  was  also 
recorded  on  the  tape.  This  permitted  precise  synchronization  of  the 
tapes  from  the  aircraft  and  from  Los  Banos,  and  prevented  dropouts  in 
the  50  kHz  signal  from  degrading  long  term  accuracy.  Voice  countdowns 
on  the  tape  allowed  exact  matching  of  the  time  between  Los  Banos  and  the 
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or  of  aircraft.  Figure  48  Inter 
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aircraft.  Tests  of  this  method  of  recording  and  the  associated  data- 


reduction  process  showed  that  a  phase  accuracy  of  better  than  ±5  deg  was 
maintained  for  the  data  channels  (where  only  a  single  frequency  component 
was  present)  and  ±10  deg  for  the  tracking  channel  (two  frequencies,  dif¬ 
fering  in  strength  by  up  to  30  dB,  were  used  for  the  tracking  system  as 
described  below).  Figure  49  is  a  block  diagram  of  the  equipment  on  the 


a ircraf t . 


The  receiving  equipment  at  Los  Banos  used  basically  the  same  tech¬ 


nique.  A  single  1/4  wavelength  vertical  antenna  was  used.  Again,  all 


conversion  frequencies  were  derived  from  a  cesium  beam  frequency  stand¬ 


ard  .  The  same  data  recording  scheme  was  used . 


CW  transmissions  from  Bearden,  Arkansas  were  provided  by  a  crystal 


frequency  standard,  a  frequency  synthesizer  and  a  1000  watt  linear  ampli¬ 


fier.  The  antenna  used  was  a  vertical  monopole.  The  frequency  stability 


of  the  crystal  frequency  standard  is  specified  as  10  for  averaging 


times  from  1  sec  to  1  hr.  This  can  introduce  an  uncorrected  doppler  fre¬ 


quency  of  up  to  0.0025  Hz  for  a  typical  pass.  The  peak -measured  doppler 
(dP/dt)  for  all  flights  ranged  from  0.05  Hz  to  0.2  Hz.  Errors  in  rela¬ 


tive  doppler  between  Los  Banos  and  the  aircraft  are  determined  by  the 


stability  of  the  cesium  standards  at  Los  Banos  and  on  the  aircraft.  These 


errors  were  of  the  order  of  0.00025  Hz. 


C.  AIRCRAFT  TRACKING  SYSTEM 


An  accurate  determination  of  wave-front  shape  from  a  moving  aircraft 


requires  an  accurate  means  of  determining  the  location  of  the  aircraft, 


The  accuracy  with  which  the  position  of  the  aircraft  can  be  determined 


sets  a  limit  on  the  accuracy  of  the  measurements  which  can  be  made.  A 


tracking  system  was  used  for  these  experiments  which  provided  continuous 


measurements  of  the  position  of  the  aircraft.  The  estimated  accuracies 


for  this  tracking  system  are  discussed  below  and  summarized  in  Table  1. 
The  tracking  system  consisted  of  a  CW  doppler  system,  an  aneroid 


differential  altimeter,  and  an  aerial  camera.  The  CW  doppler  system  mea¬ 


sured  changes  in  distance  to  the  aircraft  from  two  ground  points  within 
a  direct  line-of -sight  of  the  aircraft.  Two  doppler  transmitting  sites 


SEL-70-037 


Figure  49.  Block  diagram  of  major  equipment  of  the 
a ircraf  t . 
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were  used  (see  Fig.  46),  one  located  near  the  southern  end  of  the  flight 
path  (Los  Banos),  and  one  to  the  east  of  the  flight  path  in  the  Sierra 
foothills  (Bear  Mountain).  Both  transmission  paths  to  the  aircraft  were 


unobstructed  along  the  entire  flight  path.  Frequencies  near  25  MHz  were 


used  for  reasons  of  simplicity  and  economy  in  the  transmitting  and  re¬ 


ceiving  equipment.  The  frequencies  of  the  two  transmitters  were  offset 


by  3.0  kHz.  The  signals  were  received  in  the  aircraft  using  the  same 


vertical  antenna  that  was  used  for  receiving  the  ionospherically-propa- 


gated  signals.  Both  tracking  signals  were  received  by  a  single  receiver 


having  a  6  kHz  bandwidth.  The  audio  output  of  this  receiver  (at  3  and  6 
kHz)  v/as  recorded  on  one  track  of  the  analog  tape  recorder.  All  trans¬ 


mitting  and  receiving  mixing  frequencies  for  the  tracking  system  were 
synthesized  from  a  single  frequency  standard  at  each  location.  The  fre¬ 


quency  standards  in  the  aircraft  and  at  Los  Banos  were  cesium-beam  stan¬ 


dards.  A  crystal  frequency  standard  was  used  at  the  Bear  Mountain  trans¬ 


lating  site,  in  order  to  remove  the  effect  of  frequency  fluctuations 


in  the  crystal  frequency  standard,  the  signal  transmitted  from  Bear  Moun¬ 


tain  was  also  received  at  Los  Banos,  using  a  receiver  having  all  conver¬ 


sions  synthesized  from  the  Los  Banos  cesium  standard,  and  recorded  for 


use  during  the  tracking  data  reduction, 


The  specified  accuracy  and  long-term  stability  of  the  cesium-beam 
frequency  standards  is  ±1  X  10  11 .  Short-term  stability  is  such  that 


the  accumulated  phase  error  never  exceeds  that  due  to  the  accuracy  spe¬ 


cification  by  more  than  1  milliradian  at  25  MHz.  Measurements  made  on 


the  two  cesium  standards  that  were  used  indicate  a  relative  drift  of 


less  than  1  |asec/day  between  them.  This  is  within  the  specification 


value.  A  typical  pass  over  the  flight  path  requires  about  1/2  hr.  Dur¬ 


ing  this  period,  the  total  phase  change  is  approximately  25  X  10  X 


2 it  x  1800  rad.  The  phase  error  expected  would  then  be 


6 

25  X  10  X  2tc  x  1800  X  10  =  2.8  rad 


For  a  frequency  of  25  MHz  (A  =  12)  this  corresponds  to  a  distance  error 


of  5  m  over  a  30  min  flight,  and  proportionally  less  over  shorter  portions 


of  a  flight 
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Changes  in  the  height  of  the  aircraft  were  measured  with  a  precision 
aneroid  altimeter.  It  was  used  in  a  differential  mode  with  a  full-scale 
output  of  ±40  m.  In  this  mode  the  primary  source  of  instrument  error  is 
hysteresis,  which  is  specified  at  ±3  m.  Corrections  for  variations  in 
air  pressure  with  time  and  along  the  flight  path  were  made  when  necessary, 
using  temperature-corrected  barometric  readings  from  the  Sacramento  and 
Stockton,  California,  Weather  Bureau  stations.  These  readings  were  taken 
at  half  hour  intervals  during  the  period  of  the  flights.  For  a  majority 
of  the  flights  no  correction  was  made,  as  the  indicated  correction  was 
less  than  10  m.  Because  the  elevation  angle-of -arrival  for  the  ionospher- 
ically  propagated  signals  from  Bearden  is  low  (~  10  deg),  errors  in  alti¬ 
tude  affect  the  phase-path  measurements  by  less  than  1/5  as  much  as  do 
errors  in  horizontal  position. 

The  three  measurements  made  by  the  doppler  system  and  the  altimeter 
are  sufficient  to  determine  the  location  of  the  aircraft,  continuously 
if  an  initial  location  can  be  determined  from  other  means.  For  the  mea¬ 
surements  involved  in  these  experiments,  although  it  is  important  to  know 
accurately  only  changes  in  the  great  circle  distance  to  the  transmitter 
and  changes  in  altitude,  the  geometry  of  the  doppler  system  requires  that 
a  fixed  point  be  used  as  a  starting  reference.  The  error  in  change  of 
distance  will  be  about  half  the  error  in  the  location  of  the  fixed  point 
for  the  transmitter  locations  employed.  (If  only  portions  of  the  flight 
path  near  the  center  are  used,  where  the  angle  between  the  flight  path 
and  the  transmission  from  Bear  Mountain  is  reasonably  large,  less  error 
in  change  of  distance  will  result.) 

Photogramme trie  reduction  of  the  aerial  photographs  taken  over  two 
"photo  points"  (Fig.  46)  on  the  flight  path  (a  stereo  pair  at  each  point) 
was  used  to  establish  the  fixed  point  required  by  doppler/altimeter  sys¬ 
tem,  and  to  check  the  closure  error  of  the  tracking  system.  The  stated 
accuracy  of  the  photographic  position  determination  is  ±50  m  in  the  hor¬ 
izontal  plane,  and  ±100  m  in  altitude. 

The  final  accuracy  of  the  tracking  system  is  estimated  in  terms  of 
the  effects  of  various  errors  on  relative  position,  rather  than  absolute 
position  determination  which  is  not  important  for  the  experiments.  Table 
1  summarizes  the  various  sources  of  error  in  the  tracking  system,  and 
their  contribution  to  relative  position  error. 
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TABLE  1 


AIRCRAFT  TRACKING  SYSTEM  ERRORS 
(30  Min.  Flight  Time) 


Estimated  ms 


Error  Source  Error  (m) 


Frequency  Standards  5.0 
Altimeter  (3.0/5)  0.6 
Air  Pressure  Changes  (10.0/5)  2.0 
Photogrammetry — Horizontal  (50/2)  25.0 
Photogrammetry — Vertical  (100/10)  10.0 
Tracking  Frequency  Ground  Reflections  (A/6)  2.0 
Phase  Measuring  System  (Insignificant)  — 

TOTAL  27.4 


As  a  check  on  long-term  error,  photographic  position  determinations 
were  made  at  both  ends  of  the  flight  path  for  four  flights.  These  were 
compared  with  the  positions  determined  by  the  tracking  system.  Closure 
errors  were  20,  25,  35,  and  60  m  for  these  flights.  No  systematic  er¬ 
rors  were  visible  on  the  data.  It  is  therefore  believed  that  tracking- 
system  errors  contributed  an  rms  error  of  about  10  per  cent  to  the  ex¬ 
perimental  data  from  the  aircraft. 


D.  DATA  PROCESSING 

Each  pass  over  the  flight  path  resulted  in  one  analog  tape  each  from 
the  aircraft  and  from  Los  Banos.  Table  2  lists  the  channel  assignments 
for  these  tapes.  The  first  step  in  the  data  reduction  was  to  play  back 
the  tapes,  and,  using  a  Scientific  Data  Systems  Sigma  V  computer  with  an- 
alog-to-digital  converter,  convert  the  audio  signals  to  a  digital  tape 
which  represented  the  total  number  of  cycles  (including  fractions)  from 
the  start  of  the  tape.  This  was  done  for  each  of  the  six  audio  signals 
on  the  tapes.  Sampling  of  the  signals  was  clocked  by  the  50  kHz  or  25 
''Hz  tones  on  the  tapes.  The  one  pps  signal  was  also  sampled  to  indicate 

the  start  of  each  second,  in  the  event  that  dropouts  occurred  in  the  50 

kHz  signal.  A  test  was  made  during  sampling  to  assure  that  between  0.8 

and  1.2  seconds  of  real  time  elapsed  between  occurrences  of  the  1  pps 
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TABLE  2 


DATA  RECORDED  ON  ANALOG  TAPES 


Aircraft  Tapes 

—  30  in. /sec  —  Direct  record  (except  FM  where  marked 

with  asterisk) 

Track 

1. 

Tracking  Receiver  Output  -  3  Sc  6  kHz 

2.* 

Voice  commentary  Sc  countdowns,  Sc  camera  exposure  tone 

3. 

1  pps  from  cesium  standard 

4.* 

Differential  altimeter  - ±  40  m 

5  . 

50  kHz  sine  wave  from  cesium  standard  and  divider 

6. 

Horizontal  antenna  data-receiver  output  -  3  kHz 

7  . 

Vertical  antenna  data-receiver  output  -  3  kHz 

Los  Banos  Tapes^ —  15  in  ./sec 

Track 

1 

Beam  Mountain  tracking  correction  signal  -  3  kHz 

2.* 

Voice  commentary  and  countdowns 

3. 

1  pps  from  cesium  standard 

5. 

25  kHz  sine  wave  from  cesium  standard  and  divider 

7. 

Data  receiver  output  -  3  kHz 

4  &  6. 

Not  used . 

signal  from  the  tape.  This  insured  that  seconds  would  not  be  "lost"  if 
there  were  dropouts  in  the  1  pps  channel.  Finally,  the  altimeter  channel 
and  its  calibration  were  sampled  each  second. 

The  resulting  digital  tapes  were  listed,  and  a  test  was  made  to  de¬ 
termine  that  the  doppler  frequency  each  second  was  within  1/2  Hz  of  the 
doppler  frequency  for  the  preceding  second.  If  not,  the  data  were  flag¬ 
ged  at  that  point  to  indicate  a  tape  or  signal  dropout.  Between  0  and 
40  points  for  each  tape  (in  a  total  of  about  1500  points)  were  flagged. 

An  editing  program  was  then  used  to  correct  flagged  points,  align  all 
tape  playbacks  to  a  common  starting  point,  and  convert  from  cycles  to 
distance.  The  latter  conversion  required  inversion  of  some  of  the 
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frequencies  (where  the  receiver  used  an  odd  number' of  inverting  conver¬ 
sions),  and  multiplication  by  the  particular  wavelength  for  the  frequency 
involved.  The  resulting  tape  was  then  processed  to  obtain  the  required 
results  in  terms  of  changes  in  phase  path  length  '^.o  the  transmitter. 

Three  of  the  six  audio  frequencies  on  the  analog  tapes  result  in  two 
records,  describing  changes  in  distance  between  the  aircraft  and  Los 
Banos  and  between  the  aircraft  and  Bear  Mountain.  The  altimeter  samples 
result  in  a  record  of  changes  in  altitude  of  the  aircraft.  (Corrections 
for  changes  in  barometric  pressure  are  also  included  in  the  editing  pro¬ 
gram.)  These  three  records  are  used,  together  with  an  initial  point  de¬ 
termined  from  photography,  to  calculate  the  location  of  the  aircraft  each 
second,  and  to  calculate  the  great  circle  distance  to  the  transmitting 
site  each  second. 

The  remaining  three  audio  frequencies  produce  records  representing 

the  change  in  phase  path  length  of  the  ionospherically  propagated  signal 

✓ 

received  at  Los  Banos  and  on  the  aircraft  horizontal  and  vertical  anten¬ 
nas.  The  data  from  the  aircraft  are  then  corrected  for  the  effect  of 
changes  in  aircraft  distance  from  the  transmitter  and  changes  in  altitude. 
The  resulting  records  for  the  aircraft  represent  measurements  comparable 
to  the  Los  Banos  measurements.  Finally  a  constant  doppler  (linear  phase- 
path  change)  term  is  subtracted  from  all  the  data  for  a  given  flight. 

-The  quantity  subtracted  is  the  average  doppler  at  Los  Banos  for  a  period 
of  a  half  hour  to  an  hour,  including  the  time  of  the  flight.  The  result¬ 
ing  records  show  changes  in  phase-path  length  between  the  transmitter  and 
Los  Banos,  and  between  the  transmitter  and  the  aircraft,  with  a  linear 
term  removed  from  both.  Typical  records  are  included  in  Chapter  II.  The 
difference  in  doppler  at  varying  separations  of  receiving  sites  is  deter¬ 
mined  by  measuring  the  slopes  of  the  two  curves  at  various  times  and  cor¬ 
responding  separations. 


E.  LATER  LOS  BANOS  MEASUREMENTS 


The  measurements  made  during  the  winter  and  spring  of  1970  involved 
only  the  Los  Banos  receiving  site  and  the  Bearden  transmitting  site. 

The  basic  measurement  was  the  same,  but  the  equipment  used  was  somewhat 

129  SEL-70-037 


different.  The  transmitting  rfite  transmitted  two  synthesized  frequencies 
The  receiving  site  was  equipped  with  a  small  computer  with  an  analog-to- 
digital  converter  and  a  digital  tape  recorder.  These  were  used  to  con¬ 
siderably  simplify  the  data  processing.  Receiving  antennas  were  availa¬ 
ble  which  were  located  on  a  line  perpendicular  to  the  direction  of  the 
transmitting  site.  These  were  eight  vertical  monopoles  spaced  360  m  be¬ 
tween  each  element.  A  maximum  of  three  elements  was  used  at  any  one  time 
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APPENDIX  B 


THE  EFFECT  OF  CHANGES  IN  RAY  REFLECTION  HEIGHT 


This  appendix  discusses  the  possibility  that  the  relation  which  was 
derived  between  the  measured  changes  in  phase-path  length,  and  the  change 
in  refractive  index  and  electron  density  which  are  required  to  produce 
them,  may  not  be  accurate.  Two  distinct  simple  models  of  the  ionosphere 
are  commonly  used  in  considering  the  source  of  doppler  effects  on  iono- 
lspherically  propagated  waves  (Davies  ’65).  The  first  assumes  that  vari¬ 
ations  in  refractive  index  occur  along  straight  (non-deviating)  portions 
of  the  ray  paths  (i.e.,  not  near  the  reflection  height).  This  is  the 
assumption  made  in  the  two-dimensional  model  used  in  Chapter  III.  The 
second  commonly-used  model  assumes  that  if  there  are  ionospheric  varia¬ 
tions  near  the  reflection  height  (where  in  fact  it  appears  that  the  ob¬ 
served  irregularities  are  located),  these  can  be  analyzed  in  terms  of  a 
change  in  height  of  a  reflecting  mirror.  This  latter  idea  implies  that 
an  irregularity  may  change  the  height  at  which  a  ray  is  reflected.  The 
effect  of  such  a  change  in  reflection  height  is  now  to  be  considered. 

It  might  be  anticipated  that  when  an  irregularity  disturbs  the  path 
of  a  ray  in  the  ionosphere,  causing  a  change  in  the  reflection  height  of 
the  ray,  the  measured  change  in  path  length  may  be  greater  than  that 
which  would  occur  for  a  ray  passing  directly  through  the  irregularity 
when  the  latter  is  so  located  as  not  to  affect  the  reflection  height — 
i.e.,  when  the  irregularity  is  on  the  ray  path  but  not  near  the  region 
of  significant  ray  bending.  Figure  50  shows  the  effect  of  a  change  in 
reflection  height  caused  by  the  presence  of  an  irregularity  near  the 
region  of  significant  ray  bending.  From  Eqs.  (5)  and  (4)  we  have 


dp  =  r 

dt  -4ay  path  ^ 


ds 


(B.l) 


It  will  now  be  shown  that  the  change  in  reflection  height  does  not  affect 
the  relations  which  were  used  in  determining  changes  in  refractive  index 
and  electron  density. 
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Figure  50.  Ionospheric  ray  paths  showing  a  change  in  reflection 
height  caused  by  an  irregularity  near  region  of  sig¬ 
nificant  ray  bending. 
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Consider  first  a  non-moving  irregularity  described  by 

^  =  f (t)  •  g(x, y, z) 

M 

In  this  expression  g(x,y,z)  describes  a  fixed  shape  for  the  irregular¬ 
ity,  and  f(t)  describes  some  time  function.  Let  the  irregularity  be 
isolated  from  the  slowly  varying  ionospheric  background,  i.e.,  let  it  be 
immersed  in  a  medium  of  constant  refractive  index.  At  any  time  t,  de¬ 
fine  a  geometrical  ray  path  within  the  irregularity  as  the  path 

which  exhibits  the  maximum  doppler  at  the  time.  That  is, 


(B .  2 ) 


✓ 

Now  at  any  instant  t,  the  rate  of  change  of  phase  path,  dP/dt,  along 

the  geometrical  path  G  will  correspond  to  the  maximum  possible  value 

m 

of  the  integral  in  Eq.  (B.l).  If  the  irregularity  is  now  immersed  in 
the  ionosphere,  the  only  portion  of  any  ionospheric  ray  path  which  can 
contribute  to  a  change  in  the  value  of  the  integral  is  the  portion  (if 
any)  which  lies  within  the  irregularity,  since  outside  the  irregularity, 
dp/dt  =  0.  Therefore,  at  any  instant,  the  doppler  over  the  ionospheric 
path  must  be  less  than  or  equal  to  the  doppler  along  Gm  within  the 
isolated  irregularity.  For  simple  shapes  of  irregularities  of  arbitrary 
strength  (excluding  the  possibility  of  total  ray  reflection)  and  for  ar¬ 
bitrary  shapes  of  weak  irregularities,  we  can  show  that  Gm  is  a  fixed 
path.  In  the  case  of  a  spherically  symmetric  irregularity,  G^  is  any 
diameter  of  the  irregularity.  For  an  ellipsoidal  irregularity,  G^  is 
the  major  axis  of  the  irregularity.  In  the  case  of  weak  irregularities, 
where  in  this  case  'weak'  implies  that  rays  are  not  displaced  within  the 
irregularity,  but  only  deflected  in  angle,  all  ray  paths  within  the  ir¬ 
regularity  are  straight  lines,  and  every  straight  line  is  a  ray  path. 

G  is  then  a  fixed  straight  line  defined  by 
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g(x,y,z)  =  maximum 


(B .  3  ) 


-'STRAIGHT 

LINE 

The  irregularities  which  were  described  based  on  the  ionospheric 
measurements  of  this  work  satisfy  the  'weak'  approximation  with  little 
error.  The  result  states  that  regardless  of  the  smooth  bending  geometry 
of  the  ionosphere  itself,  irregularities  produce  a  doppler  which  is  at 
every  instant  less  than  or  equal  to  the  doppler  that  would  be  measured 
for  a  path  in  the  irregularity  when  the  latter  is  isolated  from  the 

ionosphere.  Therefore,  the  total  measured  change  in  phase-path  length 
when  the  irregularity  is  in  the  ionosphere  will  be  less  than  or  equal  to 
the  change  in  phase-path  length  along  a  path  G^  through  the  isolated 
irregularity.  So,  despite  the  fact  that  an  irregularity  may  change  the 
reflection  height  of  a  ray  in  the  ionosphere,  and  therefore  also  change 
the  phase-path  length  of  the  ray  outside  the  region  of  the  irregularity, 
the  total  measured  change  will  be  no  greater  than  the  change  which  would 
be  predicted  from  propagation  through  an  isolated  irregularity,  having 
the  same  profile,  which  is  not  close  to  the  reflection  height. 

If  the  irregularity  is  assumed  to  move,  little  is  changed  in  the 
argument.  Since  the  total  ionospheric  ray  path  at  a  given  time  is  not 
a  function  of  its  history,  the  integral  of  dP/dt  will  be  the  same 
whether  the  irregularity  reached  its  present  position  by  moving  there 
or  by  developing  there.  This  leads  to  the  conclusion  that  the  method 
used  in  relating  AP  to  A|j  is  valid. 
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The  primary  purpose  of  this  investigation  was  to  determine  the  nature  and 
characteristics  of  a  class  of  large,  weak  ionospheric  irregularities  which  in  cer¬ 
tain  circumstances  significantly  affect  long-distance  HF  radio-wave  propagation. 

A  further  purpose  was  to  determine  the  relative  practical  importance  of  the  dif 
fraction  produced  by  these  inhomogeneities,  in  relation  to  the  influence  of  other 
mechanisms  such  as  Faraday  rotation  and  scattering  by  small  irregularities. 

Measurements  of  quasi-period ic  variations  in  phase-path  length  have  been  made 
on  a  2600  km  one-hop  lower-ray  F-region  path.  Variations  of  8  to  75  wavelengths 
are  regularly  observed  during  winter  daylight  hours.  They  have  a  quasi  period  of 
about  20  minutes.  Correlations  of  measurements  made  at  spaced  receiver  locations 
indicate  that  these  variations  in  phase-path  length  are  produced  by  irregularities 
having  a  horizontal  scale  size  of  about  35  km.  Multiple  frequency  measurements 
show  a  vertical  scale  size  of  about  15  km,  with  a  maximum  in  the  lower  F  region, 

A  lower  limit  of  65  m/sec  is  determined  for  the  horizontal  velocity  of  the  irreg¬ 
ularities.  A  downward  vertical  velocity  of  30  to  40  m/sec  is  observed.  The  irreg 
uiarities  limit  achievable  bandwidth  and  spatial  resolution  for  one-hop  HF  iop.o 
spherically  propagated  waves  under  quiet  ionospheric  conditions. 

^  ”  a  model  of  the  irregularities  has  been  devised  and  has  been  used  to  predict 

the  signal-strength  fluctuations  which  would  result  from  propagation  through  a 
medium  described  by  the  model.  Two  classes  of  signal-strength  behavior  are  pre¬ 
dicted  and  observed  experimentally,  depending  on  the  relative  magnitude  of  the 
phase-path  variations  within  the  limits  mentioned  above.  Signal-strength 
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